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EXPERIMENTAL GENERATION OF PLASMA ALFVEN WAVES* 


T. K. Allen,t William R. Baker, Robert V. Pyle, and John M. Wilcox 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received April 2, 1959) 


Alfvén’ has postulated the existence of hydro- 
magnetic waves to explain certain sunspot phe - 
nomena. Some rather complex experimental ob- 
servations have been reported.? We have at- 
tempted to generate such a wave at one end of a 
cylindrical plasma, and to detect the transmitted 
wave at the other end. This hydromagnetic wave 
guide consists of a conducting cylinder, 34 in. 
long and 53 in. in diameter, filled with hydrogen 
gas to a pressure of 100 microns, and mounted 
ina uniform axial magnetic field. At each end 
of the cylinder there is an insulator and a small 


coaxial electrode, 2 in. long and 2 in. in diameter. 


The velocity of propagation of these waves can be 
derived® from the formula v 4 = c/ VK, where K, 
the dielectric constant of a plasma, is equal to 
1+(4npc*/B*); p is the plasma density, B the 
magnetic field, and c the velocity of light. Since 
K>>1, we have 


v4 7 B/(4p). (1) 


The tube is preionized by discharge of a “slow” 
condenser bank (45 yf) between the two central 
electrodes. This current reaches a maximum 
value of 50 kiloamperes in 17 ysec, and then 
“tings” a few times. A fast condenser bank 
(l.2 uf) is then discharged between one center 
electrode and the outer cylinder. This current 
Teaches a maximum value of 25 kiloamperes in 
0.7 usec and also rings a few times. The radial 
electric field thus produced and the steady axial 
magnetic field generate a torsional hydromag- 


netic wave. When this wave arrives at the far 
end of the cylinder, the radial electric field E, 
=- Vo xB, produces a voltage difference between 
the center electrode and the outer cylinder. In 
Fig. 1A, the bottom trace of the dual-beam os- 








FIG. 1. Dual-beam oscilloscope traces of the trans~ 
mitted wave form (top) and the driving wave form 
(bottom). The time scale is 1sec/cm. The top volt- 
age scale is 500 volts/cm, and the bottom is 2500 v/cm. 
The hydrogen gas pressure is 100 microns. The mag- 
netic fields are 9.0, 8.1, 7.2, 6.3, 5.4, and 4.5 
kilogauss in pictures A through F, respectively. Pic- 
ture G shows the “prompt” transmission at low ion 
density. 
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cilloscope shows the driving wave form, while 
the top trace is the transmitted wave form. The 
time scale is 1 usec/cm. The transmission 

time of the wave is measured from the first peak 
of the driving wave form to the first peak of the 
transmitted wave form. From this transmission 
time and the known length of the cylinder we have 
a direct measurement of the wave velocity. 

The operation of the measuring circuits can be 
checked by discharging the fast condenser bank 
when there is very slight ionization in the tube. 
In this case the transmitted wave arrives at the 
far end, as would be expected, with no meas- 
urable delay (Fig. 1G). 

The change in the transmitted wave form as 
the magnetic field is varied is shown in Fig. 1, 
A through F. The transmission improves with 
the larger magnetic fields. The experimental 
velocities are plotted versus magnetic field in 
Fig. 2, where the solid curve represents the 
theoretical value given by Eq. (1). One might 
expect the experimental velocities to be some- 
what larger than the theoretical velocities, as 
follows. From the nature of the preionizing 
pulse we would expect the ionization to decrease 
at large radii. The wave velocity would then be 
a function of the radius. This would increase the 
average velocity, and also make the period of 
the transmitted waye form greater than that of 
the driving wave form, as observed in Fig. 1. 

The measured velocity decreases as the gas 
pressure is increased, but the variation is not 
as large as predicted by Eq. (1). Since the en- 
ergy from the preionizing condenser bank is ap- 
proximately constant, the percent ionization de- 
creases at higher pressures and causes a devia- 
tion from the theoretical velocity. If the ioniza- 
tion parameters are not optimized, the measured 
velocity is larger than the Alfvén velocity, and 
in the limit of very low ionization the “prompt” 
transmitted pulse of Fig. 1G is observed. 

At 12.6 kilogauss the peak of the driving wave 
form is about 1200 volts and the peak of the 
transmitted wave form is about 500 volts. This 
factor of 2.4 in the attenuation could result from 
several causes. The radially dependent wave 
velocity mentioned above can significantly de- 
crease the transmitted wave peak. Sheath effects 
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FIG. 2. Wave velocity vs magnetic field. The circles 
are experimental points and the solid curve is the 
Alfvén velocity. The hydrogen gas pressure is 100 
microns. 


may cause only part of the driving voltage to be 
effective in actually driving the wave. Some of 
the wave energy may be lost to the insulators. 
Finally there is the attenuation of the wave pass- 
ing through the plasma. 

This work was initiated after discussions with 
Harold Furth. We wish to thank Richard O’Sulli- 
van, Pierre Pellissier, Louis Biagi, and Gerald 
Wilson for able assistance. 





"Work done under the auspices of the U. S. Atomic 
Energy Commission. 

Wisitor from the Atomic Energy Research Establish- 
ment, Harwell, England. 

‘H. Alfvén, Arkiv Mat. Astron. Fysik 29B, No, 2 
(1943). 

"See T. G. Cowling, Magnetohydrodynamics (Inter- 
science Publishers, Inc., New York, 1957), Chap. 3, 
for a review to 1957; also Sawyer, Scott, and Stratton, 
Phys. Fluids 2, 47 (1959). 

3... Spitzer, Jr., Physics of Fully Ionized Gases 
(Interscience Publishers, Inc., New York, 1956). 
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CONNECTION BETWEEN TRUE EFFECTIVE MASS AND OPTICAL 
ABSORPTION IN INSULATORS 
Vinay Ambegaokar and Walter Kohn 
Carnegie Institute of Technology, Pittsburgh, Pennsylvania 
(Received April 10, 1959) 
In the independent-electron model of an insula - pansion 
i the effective mass 
tor with cubic symmetry, the effective m E,(® =E,(0)+ (2 /2m*) i+... (6) 


associated with a simple band n =0 is given by 
the well-known sum rule’ 


m 
m* =1 ~Zigton ’ (1) 


where the oscillator strength f,,, is defined by 


2\p |" 
fon =m BOE,” 
n ° 


(2) 


Consider, in the same model, a small number, 
tm, of electrons per unit volume in the band 0. 
Their contribution to the real part of the com- 
plex conductivity is 


w 2 
du)= 2 fon MOEy ~ E, - He) 


+ (Ey, -E,+hw)], (3) 
where 
wy =4ange"/m. (4) 
Thus one has the following connection: 


m 8 ¥ 
— Se 1- orf ow) dw. (5) 


We have studied a more realistic model of an 
insulator in which the Coulomb interactions be- 
tween all electrons are included to all orders of 
perturbation theory,” and have found that a rela- 
tion very similar to (5) still holds. 

In this model the effective mass of a free car- 
rier is defined as follows. The many-electron 
wave functions of the insulator with one free 
carrier can be characterized by a wave vector K 
and other quantum numbers n. The energy spec- 
trum of this system will be denoted by £,,(K). 

Let E,(0) be the lowest energy. Then the many- 
article effective mass m* is defined by the ex- 


Now let o,(w) and o(w) be the real parts of the 
conductivity of the perfect insulator, without any 
free carriers, and of the insulator, with n, car- 
riers per unit volume in the lowest “conduction 
band, ” respectively. Then one finds the follow- 
ing exact® relationship: 


Ly zr f oe (7) 


Under favorable circumstances, when m* is 
very small, this relation may enable one to de- 
termine the effective mass by optical measure- 
ments. The main effect of the presence of car- 
riers in such a case is to reduce the optical ab- 
sorption of the perfect insulator near the absorp- 
tion threshold. The amount of this reduction 
gives an estimate of m* by means of Eq. (7).* 

We should like to take this opportunity to thank 
Dr. M. Suffczynski for many stimulating dis- 
cussions. 





‘See, e.g., N. F. Mott and H. Jones, Theory of the 
Properties of Metals and Alloys (Oxford University 
Press, Oxford, 1936), Chap. II. 

*w. Kohn, Phys. Rev. 105, 509 (1957); 110, 857 
(1958); Vinay Ambegaokar and Walter Kohn, Bull. Am. 
Phys. Soc. Ser. Il, 4, 276 (1959). It should be re- 
marked that in all this work the lattice is taken to be 
rigid. 

‘In particular, all exchange effects are completely 
contained in this relation. This may be contrasted 
with the analyses of E. O. Kane, J. Phys. Chem. 
Solids 6, 236 (1958), and J. C. Phillips, J. Phys. 
Chem. Solids 7, 52 (1958). 

‘Compare H. Y. Fan and G. W. Gobeli, Bull. Am. 
Phys. Soc. Ser. II, 3, 111 (1956); W. G. Spitzer and 
H. Y. Fan, Phys. Rev. 106, 882 (1957). 
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NUCLEAR MAGNETIC RESONANCE IN SUPERCONDUCTING TIN* 
G. M. Androes and W. D. Knight 
University of California, Berkeley, California 
(Received April 3, 1959) 


The difficulties of interpretation that plagued 
previous nuclear magnetic resonance experiments 
in superconducting metals’’? arose mainly be- 
cause the samples contained widely varying par- 
ticle sizes. Thus, a significant fraction of the 
sample remained in the normal state at tempera- 
tures < T, in the fields employed. Since the two 
different types of particles contributed resonance 
lines which overlapped, it was difficult to sepa- 
rate the part arising from the absorption in the 
superconducting material. In order to avoid this 
difficulty, a sample has been built up by evapo- 
rating, alternately, layers of tin and a dielectric 
(Nylon). The metal in the sample (approximately 
one gram) has the £ tin structure, and is uni- 
formly divided into small platelets with diameters, 
in the plane of the evaporated layers, of ~140 A, 
and thickness ~40 A. The critical temperature 
is found to be 3.712 +0.010°K, and the critical 
field at T=0 is 25+3 kilogauss. Experiments 
have been performed with the sample in the super- 
conducting state in fields up to 8.8 kilogauss. 

The resonance has an almost Gaussian shape, 

a width which varies linearly with the magnetic 
field, and, aside from the usual T~' variation of 
intensity, no observable temperature dependence. 
That is, to the accuracy of measurement (signal 
to noise ratio of <10 at 4.2°K), the line shape and 
width are the same above, in, and below the su- 
perconducting transition. Work on other tin sam- 
ples containing somewhat larger particles indi- 
cates that the line width is also inversely propor- 
tional to particle size, and in the final sample it 
is ~4.5 times the width for “bulk” tin (taken as 

Vy -V My in the same field. These characteristics 
may be explained by considering, in conjunction 
with the theory of Bloembergen and Rowland, * 
that the conduction electrons near the surface of 
a piece of metal do not see the regular tin lattice 
symmetry. In the platelets in the sample, half 
or more of the atoms are within 10 A of the sur- 
face. The variety of the surface positions in re- 
lation to the crystal axes, and the effects of 
some sort of random particle-size distribution 
account for the shape of the line. 

The nuclear magnetic resonance shift in the 
evaporated metal at 4.2°K was found to be the 
same as the isotropic shift in “bulk” tin. But 
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comparison with several diamagnetic salts‘ 
yields values of the shift in normal metal Tanging 
from 0.63 to 0.81%, depending on the salt. There. 
fore, a quantity D=(v mm, v,)/ v,, Will be used to 
describe the behavior of the resonant frequency, 
Ysos in the superconducting state. Here Vy is 
the resonant frequency in the normal metal at 
4.2°K. Figure 1 is a plot of D vs the reduced 
temperature, ¢(H)=7/T,(H), for an applied field 
of 1.2 kilogauss. The curve extrapolates to 
D,~0.2% at T=0. Whether or not this extrapol- 
tion is valid, it is clear that the curve has not 
reached an asymptotic value at ¢(H) =0.4. 
Similar curves of D vs ¢(H) are obtained in 
fields up to 8.8 kilogauss. At ¢(H) =} they fall 
within a range of 0.02% above the 1.2 kilogauss 
value in Fig. 1. The effects of the magnetic field 
are thus small. Details will be reported later. 
We conclude that the nuclear magnetic reso- 
nance shift in superconducting tin does not ap- 
proach zero as T goes to zero. If SnCl, is used 
as the reference salt, (4H/H),._~ 9-55 %, which 

















T T T T T T T 
0 4 
o 
-0,1 5 
a 06 08 10 
t(H)=T/T(H) 
FIG. 1. Deviation of the resonant frequency in super- 


conducting tin. D= (Vso Vy)/Vy, where Voc is the 
resonant frequency in the superconducting state, and 
vm is the resonant frequency in the normal state at 
4.2°K. The experimental points are for H =1.2 kilo- 
gauss; T,(H) is 3.68°K for this field. The solid line 
is a plot of the equation ¢(H) ~cosh[a(D-d)}, a anid 
being constants. The lower curve is a plot of the cham 
in frequency of the experimental oscillator (arbitrary 
units) as the temperature was lowered through the 
transition. Any correction in D for the demagnetiz#- 
tion factor of the sample is less than 0.01%. 
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is approximately 73 % of its value in normal tin. 
This is consistent with the results of Reif’ for 
mercury. 

We wish to thank E. I. du Pont de Nemours and 
Company for providing Mylar film which we used 
infabricating the sample. P. C. Lauterbur kindly 
informed us of his results before their publica- 
tion in detail. We are grateful to Professor R. 
Wild for information about particle size obtained 
from the scattering of x-rays at small angles, 
and to Mr. G. Gordon for identifying the crystal 


structure by means of x-rays. 





"This work was partly supported by the Alfred P. 
Sloan Foundation and the U. S. Office of Naval Research. 

'F. Reif, Phys. Rev. 106, 208 (1957). 

*Knight, Androes, and Hammond, Phys. Rev. 104, 
852 (1956). 

‘N. Bloembergen and T. J. Rowland, Acta Met. 1, 
731 (1953). 

‘Pp. C. Lauterbur (private communication); P. C. 
Lauterbur and J. J. Burke, Abstracts of papers pre- 
sented at 133rd meeting of the American Chemical 
Society 1958, p. 15L. 





There are three phenomena of the same general 
character that occur in rare earth metals or so- 
lution thereof: 

(1) The reduction of superconducting transition 
temperature in dilute solutions of the rare earths 
inlanthanum,* or of rare earth compounds in 
CeRu,” or OsRu,.° 

(2) An extra resistance of the pure metals (and 
presumably of the alloys*) due to spin disorder 
scattering. 

(3) An indirect exchange coupling among the 

ion spins in analogy with the Ruderman-Kittel 
mechanism of nuclear spin-spin interaction in 
metals. This indirect exchange is thought to be 
responsible for the observed ferromagnetism 
and/or antiferromagnetism at temperatures of 

the order of 100°K.® 

All three effects are supposed to arise from the 
exchange interaction AS; on*Scond of the f-shell 


-} Spins with the spins of the conduction electrons. 


The purpose of this note is to point out the bear- 
ing which spin-orbit coupling has on the interac - 
tion inthe above three cases. Except in the case 
deuropium, discussed separately, we assume 
the multiplet splitting to be much greater than 
crystal field splitting. We consider the above 


three effects in turn. 
(t) Ina previous note by Herring’ and in a de- 


lailed paper by Suhl and Matthias,* it was shown 
that on the basis of the Bardeen, Cooper, Schrief- 





EFFECTS OF SPIN-ORBIT COUPLING IN RARE EARTH METALS, AND IN 
SOLUTIONS OF RARE EARTH METALS 


R. Brout 
Cornell University, Ithaca, New York 


and 


H. Suhl 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received April 8, 1959) 


fer theory of superconductivity, ° one expects the 
exchange interaction to lower the transition tem- 
perature of the superconducting state. This de- 
pression results from the fact that, partly because 
of the energy gap, and partly because of the spin- 
coherence properties of the BCS wave function, 
the exchange scattering lowers the free energy of 
the normal state more than that of the supercon- 
ducting state. It was shown in reference 8 that if 
spin-orbit coupling is neglected, the reduction 

is proportional to S* =S(S+ 1), while when it is 
included, the reduction is very nearly propor - 
tional tu the square of the projection of S on the J 
appropriate to the ground state. The only exception 
is europium (J=0) for which this projection van- 
ishes. In almost all other cases, the agreement 
with the observed depression of Tp was sub- 
stantially improved by using J(S-J)/* rather 

than S. The reason is as follows: The calcula- 
tion of the depression of free energy is a matter 
of second order perturbation theory. The oper- 
ator Sion‘Scond connects the ground state to 


states with one conduction electron excited and 
with J either unchanged or changed by 1. Since 


in the ground state J has its maximum or mini- 
mum permissible value, *° only the J-J-1 or the 
J-~J +1 transition need be considered. For de- 
finiteness we shall from now on consider only the 
case J~J-1. Now the state 7(J-1) is removed 
from the ground state energy by A, the multiplet 
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splitting, and 4 makes its appearance in the en- 
ergy denominator, in second order perturbation 
theory, in addition to the electron excitation en- 
ergy. Transitions within the ground state of J 
involve the latter energy only, which for the im- 
portant transition is small compared with A. 
Therefore, excitation to the state y(J-1) does 
not contribute appreciably to the energy depres- 
sion except in europium, where the total effect 
should come from this excitation. In any case, 
the fact that only J~J ~-1 transitions are involved, 
enables us to calculate both the (J|S|J) and the 
(J ISiJ- 1) matrix elements by a simple sum 
rule. We know that in any representation, in 
particular in the J representation, S* must 
equal S($+1). That is to say, 


1 JISIo) 174 1 J-11S 17) |? =S(S+1). 


But the first term is S(S+1)[cos*(S, J) | -ound state: 
hence the last term on the left-hand side must be 
S(S +1)[sin?(S, J) round state’ 

To summarize: For elements other than Eu, 
only the cos* term need be considered, since the 
energy denominator of the term involving sin’ is 
of order A~1000 cm™ whereas the important 
difference in the reductions of free energy of the 
normal and superconducting states derives from 
electron excitations of order hw Debye ~ 200 cm”, 
In the case of europium, however, only the sin* 
term survives and equals unity, so thatS may 
be replaced by the spin-only value. We have com- 
puted the corresponding depression of T,. To 
account for the data’ we require a value of A be- 
tween 0.1 and 0.2 times the Debye energy, i.e., 
20 to 40 wave numbers. Even though A for 
europium should be moderately small (perhaps 
100 to 200 wave numbers) it is not easy to see 
why it should be as small as this. Perhaps part 
of the reason is that the f shell attempts to com- 
plete itself at the expense of the outer electrons, 
leading to a larger depression in 7, due to the 
resulting larger value of the spin. 

(2) We now turn attention to spin disorder re- 
sistance. Anderson and Legvoldé have calculated 
the spin disorder resistance in rare earth metals, 
to the right of Gd inclusive by a subtraction pro- 
cedure on the raw data correcting for phonon and 
impurity scattering. For all of these elements 
A>>300°K so that real collisions involving a 
change in J are negligible (like e~4/RT). The 
only resistance-contributing matrix elements 
are those diagonal in J; hence we expect a de- 
pendence like S(S+1)cos*(S, J) =[(J+1)/J]s* for 
elements to the right of Gd. A detailed com- 
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Table I. Spin disorder resistance of rare earth me 


to the right of Gd. Column 2 gives the data referred , 
Gd, corrected by Anderson and Legvold for phonon ay; 
impurity scattering. Column 3 gives the ratio expec. 
ed on the basis of spin alone (zero spin-orbit Coupling), 
Column 4 gives the ratio expected when spin-orbit 
effects are included. 





—— 








Ratio on spin Spin-orbit 
Element “Observed” basis alone factor include 
Gd 1.00 1.00 1.00 
Tb 0.76 0.76 0.67 
Dy 0.64 0.56 0.45 
Ho 0.34 0.38 0.29 
Er 0.22 0.24 0.16 
Tm 0.09 0.13 0.08 
Yb 0.01 0.05 0.02 





parison with the estimates of Anderson and Leg- 
vold (Table I) shows that accord with experiment 
is considerably worsened by inclusion of spin- 
orbit effects. As comparison with the data on 
suppression of superconductivity was improved 
by inclusion of spin-orbit effects (and this was 
comparison with “raw” data), we feel that the 
estimate of Anderson and Legvold of the spin- 
disorder resistance is quite likely to be quanti 
tively unreliable. There is no obvious reason 
why, in one situation, the orbital moment shoul 
be quenched, and not in the other. 

We may note that Eu presents an extremely in- 
teresting case as here A.J =0 is forbidden. h- 
elastic scattering will give a spin-disorder re- 
sistance proportional to e~4/kT, and crystallix 
field effects will mix a small amount of J=2 int 
the ground state to give a small temperature- 
independent term for a certain range (no spin 
order). We know of no experimental verification 
of this prediction. 

(3) Finally we point out how the indirect spir- 
spin exchange is corrected for spin-orbit c 
We follow reference 5. Then the Fourier trans- 
form of the exchange potential is 


- 1 
F(q) =A? s(ss1)[cos*G, De eb) -€+q) 


z 1 
ahaa; a or | 
It is thus seen that the short-range part of the 
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ais the lattice spacing. Thus a change in shape 
js brought about. If ferromagnetism arises from 
delicate balance of the oscillating parts of the 
Kittel-Ruderman potential, it is quite possible 
that the above-mentioned change in shape could 
pring about qualitative changes in magnetic pro- 
perties. The statistical mechanical consequences 
are being worked out at present, both for pure 
samples and for dilute alloys. 

We are much indebted to C. Herring of the Bell 
Laboratories for his fruitful suggestions in this 


work. 





‘Matthias, Suhl, and Corenzwit, Phys. Rev. Lett. 1, 


92 (1958). 
‘Matthias, Suhl, and Corenzwit, Phys. Rev. Lett. 1, 


449 (1958). 

3Suhl, Matthias, and Corenzwit, J. Phys. Chem. 
Solids (to be published). 

‘G. S. Anderson and S. Legvold, Phys. Rev. Lett. 
1, 322 (1958). 

5M. A. Ruderman and C. Kittel, 
(1954). 

®For example: Thoburn, Legvold, and Spedding, 
Phys. Rev. 110, 1298 (1958). 

'C. Herring, Proceedings of the Kamerlingh-Onnes 
Memorial Conference on Low-Temperature Physics, 
Leiden, Holland, 1958 (unpublished). 

*H. Suhl and B. T. Matthias, Phys. Rev. (to be pub- 
lished). 

*Bardeen, Cooper, and Schrieffer, Phys. Rev. 108, 
1175 (1957). 

‘Depending on whether the element is to the right or 
left of Gd. 


Phys. Rev. 96, 99 





The wavelength of light required to excite the 
maximum photoconductivity in CdS crystals has 
been found to vary with its plane of polarization 
in respect to crystal orientation.’ Gobrecht and 


‘| Bartschat® have previously observed this effect. 


The crystals we used were obtained by vapor 
phase growth ina H,-H,S atmosphere. These 
samples were ribbon-type single crystals 5 to 
10mm long, 1 mm wide, 10 to 100 y» thick with 
the c axis in the plane of the ribbon and usually 
inclined at some angle to the long, straight edges 
df the crystal. Ohmic contacts were made to the 
crystal by evaporating indium onto its ends. 

The photoconductivity was excited by mono- 
chromatic light from a Beckman DU spectrophoto- 
meter provided with a polarizer over the exit slit 
and with the plane of the crystal perpendicular to 


pol the incident light. The photocurrent was meas- 


ured with a Keithley 200B electrometer with 45 


"} or 90 volts applied across the crystal. In these 


experiments, the photocurrent was measured as 
afunction of wavelength, with the polarizer set 
aagiven angle. The relation of the monochro- 
mator polarizer to the c crystal axis was then 
determined with a polarizing eyepiece crossed 
vith respect to polarizer. The pair of crossed 
Wlarizers was rotated to determine the extinc- 
tion angles and thus the c axis. 





Spectral response of photocurrent for these 
rystals was typical of CdS, i.e., beginning with 


POLARIZATION OF PHOTOCONDUCTIVITY EXCITATION BANDS IN CdS SINGLE CRYSTALS 


R. L. Kelly and W. J. Fredericks 
Stanford Research Institute, Menlo Park, California 
(Received April 7, 1959) 


a short-wavelength tail, rising to a peak, and 
falling off rapidly at first and then more slowly 
toward the long wavelengths.” Figure 1 gives the 
change in wavelength at the photocurrent maxi- 
mum, Ad(max), as a function of the angle of 
polarization of incident light with respect to the 
c axis for a typical crystal. In all, 27 crystals 
from 14 separately grown batches were examined 
and A\(max) was found to be 39+5 A or approxi- 
mately 0.019 ev. When the electric vector E was 
parallel to the c axis, the wavelength of exciting 
light was a minimum () ) =5069 A; and when per- 
pendicular, a maximenn a 12 =5108 A). 
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FIG. 1. Change in wavelength for photoconductivity 


maximum Ad (max) as a function of angle between inci- 
dent plane polarized light and c axis of crystal. In 

the crystal shown, A, is the wavelength of maximum 
photocurrent for excitation with light parallel to the 

C axis. The average AA=39+5 A. 
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FIG. 2. Lambe-Klick-Birman model of fluorescence 
and photoconduction in wurtzite structures. The bands 
are labeled, including crystal and spin-orbit splittings. 
The “center” state is assumed to have I, symmetry. 


This shift in the energy required for excitation 
of maximum photocurrent can be interpreted on 
the basis of Birman’s development® of the Lambe- 
Klick model.* A schematic drawing of the Lambe- 
Klick-Birman model of fluorescence and photo- 
conduction for sulfides with wurtzite structures 
is given in Fig. 2. The assumptions that A_, the 
crystal splitting, is larger than A, ,, the spin- 
orbit splitting, and that the lower I, valence 
band level could be neglected, were made by 
Birman.* . 

When the polarized incident radiation has E lic, 
(1) the crystal absorbs at vy, to produce conduc - 
tion electrons and holes in the valence band I, 
(2) the holes reach equilibrium between bands I, 
and I, (3) a bound electron and free hole re- 
combine, and (4) »,+v, is emitted. When the 
polarized incident radiation has Ec, (1) the 
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crystal absorbs at v,+ v, to produce conduction 
electrons and holes in valence band I, and I,, 
and processes (2), (3), and (4) then follow. Thys 
photoconduction bands and absorption band edge” 
are dependent on the polarization of incident ra- 
diation while the fluorescence is not, as observe; 
by Lempicki.® 

On the basis of this model, », >A}. This has 
been observed. From Dutton’s’ absorption and 
reflection measurements, the calculated separa- 
tion at 90°K is 0.015 ev. Bube,® in a study of the 
structure in the photosensitivity maximum of C4 
with nonpolarized excitation, has reported a band 
with two components separated by 0.017 ev, which 
agrees well with the 0.019 ev separation reported 
here. Since \, is determined by v,+v, and dy is 
determined by y,, then it is estimated by this 
model that A, o, >0.019 ev. (Dresselhaus has 
estimated A, , to be 0.04 ev in CdS from the 
atomic spectrum of sulfur.) 

The authors thank Stanford Research Institute 
for support of this work. 
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SUPERCONDUCTING ELECTRONIC SPECIFIC HEATS, THE “EXPONENTIAL LAW,” 
AND THE BARDEEN, COOPER, SCHRIEFFER THEORY 


H. A. Boorse 
Pupin Physics Laboratories, Columbia University, New York, New York 
(Received March 18, 1959; revised manuscript received April 13, 1959) 


The Bardeen, Cooper, Schrieffer’ (BCS) theory of the second kind, and 8 =1/kT. The above ex- 





of superconductivity evaluates the molar super- pression, valid under the condition fe,>1, is 
conducting electronic specific heat, C,., as actually the first term of an expression involving 
follows: an infinite sum of the Bessel functions noted 
3 /e.\/T.\? above for which the argument is (n8e,), n being 
Ces/T¢ =773 =) [3K, (Be,) +K;(8e,)], an integer. However, only the first term of the 
4 c. sum given by the above relation need be taken 
where yT, is the electronic specific heat in the into account. 
normal state at the transition temperature T ? If values of logio(C, s /yT ;) calculated from the 
¢, is the half-width of the energy gap in the single above equation are plotted against T/T, the re- 
particle density of states, k is the Boltzmann sult is the solid curve shown in Fig. 1. For Pe 
constant, the K’s are modified Bessel functions between about 2.5 and 6, the curve may be closely 
1.0 
g — ecs 
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FIG. 1. The reduced superconducting electronic specific heat versus T ‘¢/T for 
Al, Zn, and V, showing departure from exponential behavior at values of T,/T> 4. 
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approximated by a straight line which satisfies 
the equation C, s/T¢ =aexp(-bT,/T), where a 
and b are 8.5 and -1.44, respectively. For values 
of T,,/T>6, the curve drops below this empirical 
line. At larger values of the abscissa, between 
about 7 and 11, the curve may again be approxi- 
mated by a new straight line given by the same 
relation but with the a and 6 values altered to 26 
and -1.62, respectively. These conclusions and 
their comparison with available experimental 
measurements have been presented previously 

by the author on several occasions.” The publica- 
tion of new experimental data on Al and Zn, and 
conclusions based on these data, suggest that 
further comment may be useful. 

Figure 1 presents a graph of the BCS equation, 
showing the slight downward curvature which 
gives rise to the two “exponential” laws. The 
first of these, beginning in the T,/T interval be- 
yond about 2 or 2.5, has been verified for many 
superconductors and constitutes a strong experi- 
mental confirmation of the theory. Most of the 
available data, however, do not extend much 
beyond T,,/T = 4. It is of interest to determine 
whether or not the general behavior of the spe- 
cific heat at values of T,,/T>6 is also in accord 
with the theoretical predictions, i.e., whether it 
has a small downward curvature or some other 
behavior. In this connection, the dashed line has 
been drawn in the figure to show the T° variation 
predicted by Gorter -Casimir® two-fluid model. 
The experimental points in the graph are actual, 
as well as interpolated points taken from the 
measurements of Goodman,‘ Phillips,®”* and Zav- 
aritskii.’ 

For vanadium the C,,/yT,_, values of Goodman 
at T,/T>4 show excellent agreement with the 
BCS curve. Although supporting data are absent 
for this region, earlier measurements*® (not 
shown) for T,,/T<4 are in good agreement. Three 
sets of data are exhibited for Al. Inspection shows 
that the Zavaritskii measurements follow an ex- 
ponenttial to T/T =6. If exponential behavior is 
considered as being in agreement with theory, 
then these data emphasize again only the depar- 
ture from the law of corresponding states. On 
the other hand, the aluminum data of Goodman* 
and Phillips® show a significant upward curva- 
ture beginning at about T/T =4. Such a behavior 
suggests the existence of more states at low en- 
ergies than envisioned by the theory, or the ex- 
istence of states in the energy gap, or an asym- 
metry of the gap. 
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Another example of the upward curvature in 
the C,,/yT, curve is afforded by measurements 
on Zn, although here also the situation is some- 
what ambiguous. The data of Phillips,® confirme 
at T °/T values between 1 and 2 by Seidel and 
Keesom,” show an upward curvature starting at 
T ,/T=4. This behavior is not found in the 
Zavaritskii’ data which follow an exponential. 
The difference between these results and those 
for other superconductors seems large. This is 
also noted by Khalatnikov and Abrikosov"™ who 
present a digest of data similar to Fig. 1. A 
further uncertainty respecting Zavaritskii’s Zn 
measurements appears in an unfortunate dis- 
crepancy between the graph of his C,,, data and 
the equation given to represent them. A final 
comment on the deviation of Phillips’ data from 
an exponential seems warranted in view of the 
suggestion of Seidel and Keesom”™ that the dif- 
ference may arise from a Schottky anomaly 
brought about by the interaction of the quadrupole 
moment of the nucleus with the crystalline field 
of the lattice. In this connection it may be noted 
that the similar deviation in the two Al curves 
cannot be ascribed to quadrupole coupling since 
Al has a face-centered cubic lattice. Hence it is 
possible that in Zn only a part of the difference 
arises from this source. 

The data available on tin*»** have been analyzed 
by Biondi, Forrester, Garfunkel, and Satterth- 
waite.’* They point out that although the meas- 
urements of both Goodman and Zavaritskii indi- 
cate departures from exponential behavior, the 
uncertainties are such that it is impossible to 
specify them with suitable accuracy. For this 
reason the tin data have been omitted. 

Although a T° dependence for the superconductiy 
electronic specific heat has not been observed, 
such behavior may be inferred for thorium from 
the critical magnetic field measurements of Wol- 
cott and Hein.”* In an effort to determine whether 
or not such an extreme behavior will be observed, 
an investigation of the heat capacity is currently 
under way in this laboratory. 

To sum up, it appears that at T/T values of 
four or greater there is a definite departure of 
Ces/yT¢ from exponential behavior in several 
of the superconductors and that the extent of the 
departure varies. It thus appears that a further 
refinement of the BCS theory is necessary to 
account for the experimental results. 

I wish to thank Bruce Biavati, Herbert Leupoli, 
and Stanley Zemon for assistance with Fig. 1 and 
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for some detailed calculations. 








‘pardeen, Cooper, and Schrieffer, Phys. Rev. 108, 
1175 (1957). 

%nternational Conference on the Electronic Proper- 
ties of Metals, Geneva, New York, August 25-29, 195 
(unpublished); H. A. Boorse and A. Hirshfeld, Bull. 
Am. Phys. Soc. Ser,II, 4, 5 (1959). 

ic, J. Gorter and H. B. G. Casimir, Physik. Z. 35, 
963 (1934). 

4g, B. Goodman, Compt. rend. 244, 2899 (1957). 

‘ny, E. Phillips, Proceedings of the Fifth International 
Conference on Low-Temperature Physics and Chemis- 
try, Madison, Wisconsin, August 30, 1957, edited by 
J. R. Dillinger (University of Wisconsin Press, Madi- 
gon, 1958), and private communication. 

tn. E. Phillips, Phys. Rev. Lett. 1, 363 (1958). 

'N. V. Zavaritskii, J. Exptl. Theoret. Phys. U.S.S.R. 
34, 1116 (1958) (translation: Soviet Phys. JETP 7, 773 











(1958)). 

®Corak, Goodman, Satterthwaite, and Wexler, Phys. 
Rev. 96, 1442 (1954). 

*Worley, Zemansky, and Boorse, Phys. Rev. 99, 
447 (1955). 

°G. Seidel and P. H. Keesom, Phys. Rev. 112, 1083 
(1958). 
‘lJ, M. Khalatnikov and A. A. Abrikosov, Advances 
in Physics, edited by N. F. Mott (Taylor and Francis, 
Ltd., London, 1959), Vol. 29, p. 45. 

2G. Seidel and P. H. Keesom, Phys. Rev. Lett. 2, 
261 (1959). 

18N. V. Zavaritskii, J. Exptl. Theoret. Phys. U.S.S.R. 
33, 1085 (1957) [translation: Soviet Phys. JETP 6, 837 
(1958)]. 

“Biondi, Forrester, Garfunkel, and Satterthwaite, 
Revs. Modern Phys. 30, 1109 (1958). 

‘SN. M. Wolcott and R. A. Hein, Phil. Mag. 3, 591 
(1958). 





IMAGE OF THE FERMI SURFACE IN THE VIBRATION SPECTRUM OF A METAL* 


W. Kohn 
Department of Physics, Carnegie Institute of Technology, Pittsburgh, Pennsylvania 
(Received April 6, 1959) 


The lattice vibrations of the ions in a metal are 
partly screened by the conduction electrons. We 
shall see that this screening changes rather rap- 
idly on certain surfaces in the space of phonon 
q-vectors and that therefore on these surfaces 
the frequencies w vary abruptly with q. The cal- 
culations we have done give the result that w(@) 
is a continuous function of q but that on the sur- 
faces in question 

igrad, w(@)| =a. (1) 
The location of these surfaces is entirely deter- 
mined by the shape of the electronic Fermi sur- 
face, using a simple geometrical construction. 

To explain the physical origin of this effect let 
us first describe the conduction electrons by a 
free electron gas, with Fermi wave number kf. 
One then finds that an embedded charge distribu- 
tion, 


Pext(®) = Po Je, (2) 

induces an electronic charge density 
Pei) =-Fla)poe"™, (3) 

where 
1 k q q+2k 
F F 
rece fist? (c-fos 
"ag | *9 ake] \q-2kp| |’ 





here a, is the Bohr radius. Note that near qg=2k F 


F@)= tak (1+ gp @-2hp)a-24p) (5) 


and 


dF(q)__ 1 ; , 
4a * toad? In|q-2kpl=-@. (6) 


The last equation shows an abrupt decrease of the 
ability of the electrons to screen the embedded 
charge distribution as soon as g exceeds 2k F’ 
This is due to the fact that as long as ¢<2k,_, 
Poxt(F) causes virtual excitations of some elec- 
trons with conservation of energy while when 
q>2k F Such excitations are no longer possible 
(see Fig. 1). Now a lattice vibration of wave vec- 


[4 . | 


FIG. 1. Virtual excitations for q<2ky and q> 2p. 
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tor q produces a change of ionic charge density in § 
of the form ene 
Pion®) = D, AyexpliG+K,)-F], (7) . 
where K., are the reciprocal lattice vectors. the 
Therefore we expect an abrupt change of the re- gen 
storing force whenever q is such that, for some ala 
reciprocal lattice vector K, ; of a 
- befc 
\q+K,,| =2kp. (8) - 
On the surfaces in q-space defined by (8), one wou 
finds the singularity (1) as a consequence of (6). sing 
Next we consider noninteracting Bloch electrons Il, 
with a Fermi surface given by 
ER) ot. (9) i 2. Construction for the case of Bloch elec- | = 
Again one finds singularities of grad+w, whose 
locus is determined by the following construction do not affect the nature of the singularities, but ie 
(Fig. 2): Let I, and I, be two parallel planes in only their magnitude. It therefore appears un- . 
k-space touching the Fermi surface at K, and k,. likely that the Coulomb interactions can obliter- 
Then there exists exactly one reciprocal lattice ate the effect we have discussed, which basically 
vector K, such that the vector q, defined by reflects the sharpness of the Fermi surface. Ot 
7 = -K.k (10) The magnitude of the effect may be quite large 
sail in (very roughly of the order of percent), and its 
lies in the fundamental Brillouin zone. At this observation in lattice vibration spectra would be 
point q, Eq. (1) holds. The totality of pairs of give rather direct information about the shape of Sei 
planes [I], and Il, generate the required locus of the Fermi surface. 
singularities of the vibration spectrum by means Similar “images” of the Fermi surface may be No 
of Eq. (10). expected in spin wave spectra, when the interac- 
Finally the question arises how the Coulomb in- tion between localized spins is brought about by 
teraction between the conduction electrons affects exchange with conduction electrons. 
our conclusions. We have partly allowed for this A detailed report is in preparation. 
interaction by regarding it as a perturbation and 
summing certain important terms in the result- Not 
ing perturbation expansion of w(q). These terms * supported in part by the Office of Naval Research. 
ROLE OF SINGLE PHONON EMISSION IN LOW-FIELD BREAKDOWN | 
OF SEMICONDUCTORS AT LOW TEMPERATURES joniza 
Murray A. Lampert, Frank Herman, and M. C. Steele =a 
RCA Laboratories, Princeton, New Jersey bris 
(Received April 6, 1959) are po 
Electrical breakdown at low temperatures and cussed, and some comments are made concern- a 
at low field strengths’ has been observed in both ing breakdown data for Ge doped with deep-lying tapes 
n- and p-Ge,”?~™ but has never been seen in Si.° acceptors. Finally, several experiments are which ; 
It is the purpose of this note to correlate these proposed. for the 
observations with the known energy-band struc- The LFB observed in Ge at liquid He tempera- volved 
ture and phonon spectra of Ge and Si and there- tures is produced through impact ionization of nGe a 
by to suggest a simple, necessary condition for shallow impurities by free carriers heated by — 
LFB (low-field breakdown). Difficulties in estab- the applied field. It has been suspected by those ergies 
lishing a sufficient condition for LFB are dis- working on this problem that the absence of LFB duction 
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in Si is due to the degradation of “hot” carrier 
energy Via emission of an optical phonon. Al- 
though this explanation is certainly oversimpli- 
fied, since it does not work for n-Si, a study of 
the available data does suggest a slightly more 
general proposition: If a “hot” carrier can lose 
a large fraction of its energy through emission 

of a single phonon, either acoustical or optical, 
before it reaches the threshold energy E, for 
impurity ionization, then LFB will not occur. It 
would be natural to add: “Conversely, if no such 
single phonon is available, then LFB will occur.” 





However, the failure to observe LFB in boron- 
doped Si provides a notable, if unique, exception. 
Table I lists data relevant to the LFB problem 

for Ge and Si. In m-Si several of the energetic 
intervalley phonons which are listed have low 
enough energy to prevent LFB and we do not ex- 
pect any selection rules to cut down their number. 
On the other hand, the optical phonon energy E opt 
is too great for any long-wavelength optical 
phonon to be playing a role. Inn-Ge, presumably 
the TA (transverse-acoustical) intervalley phonon 
would prevent LFB in arsenic-doped samples 


Table I. Low-field breakdown data for Ge and Si. 











LFB Material Ey (ev)* Phonon-energy (ev)? 
Observed n-Ge Sb, P, As: Eopt + 0.037% 
0.0096 through 0.0127 
Ejy : TO=0.034, L=0.029, TA=0. 010% 
TA inoperativef 
Observed p-Ge Al, B, Ga, In: Eopt : 0.037% 
0.0102 through 0.0112 
See text p-Ge Zn: 0.03 
See text p-Ge Cu: 0.04; Cd: 0.05; Egpt : 0.037%" 
Mn, Ni, Co, Fe>0.1 
Not observed = n-Si Li, Sb, P, As: Eopt : 0.0638 
0.033 through 0.049 
Eiy : TA=0.02, ey fl 
um LO=0.05, TO=0. 06 
Evy: TA=0.012, gong I 
LO=0.061, TO=0.061 
Not observed  p-Si B: 0.045 (0.046); Eopt : 0.0638" 


Al: 0.057 (0.067); 
Ga: 0.065 (0.071); 
In: 0.16 





"Ey is the thermal ionization energy, except for the 
three numbers in parentheses (p-Si), which are optical 
ionization energies. Data are taken from W. Kohn, in 
Solid State Physics (Academic Press, New York, 1957), 





Vol. 5, p. 258. 


List of “substantial” electron energy losses that 
are possible in a scattering with emission of a single 
phonon with conservation of energy and crystal momen- 
tum, Font is the energy of the k=(0,0,0) and therefore 
(approximately) of any long-wavelength optical phonon 
which may be involved in a “nearly vertical” transition 
for the electron. Ej, is the energy of a phonon in- 
volved in an intervalley transition of the electron in 
n-Ge and m-Si. Inn-Si, a transition between any two 
valleys lying along the same axis in k-space involves 
“umklapp process. The corresponding phonon en- 
ergies are labelled Ejy"™. The minima in the con- 
duction band are taken at k= 2na~'(0.85, 0,0): G. Feher 


(private communication). TO, TA, LO, LA denote 
transverse optical, transverse acoustical, longitudinal 
optical, and longitudinal acoustical phonons, respec- 
tively. L denotes a concurrence of the LO and LA 
branches of the phonon spectrum. 

°Pelah, Eisenhauer, Hughes, and Palevsky, Phys. 
Rev. 108, 1091 (1957). 

45.°N. Brockhouse and P. K. Iyengar, Phys. Rev. 
108, 894 (1957); 111, 747 (1958); B. N. Brockhouse, 
J. Phys. Chem. Solids 8, 400 (1959). 

°Ghose, Palevsky, Hughes, Pelah, and Eisenhauer, 
Phys. Rev. 113, 49 (1959). 

C. Herring, Bell System Tech. J. 24, 237 (1955); 
C. Herring and E. Vogt, Phys. Rev. 101, 944 (1956). 

&B. N. Brockhouse, Phys. Rev. Lett. 2, 256 (1959). 

Dpalevsky, Hughes, Kley, and Tunkelo, Phys. Rev. 
Lett. 2, 258 (1959). 
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were it not for the fact that the corresponding 
transition is forbidden by a selection rule. For 
p-type Si the possibility of optical phonon emis- 
sion satisfactorily explains the absence of LFB 
for In, Ga, and Al dopings. Note that it is the 
optical, rather than the thermal, ionization en- 
ergy which is significant for the breakdown, 
since impact ionization takes place before the 
lattice can substantially relax. This remark has 
particular pertinence for the case of Al-doped 
Si. 

It is puzzling that LFB occurs in B-doped Ge 
and not in B-doped Si. From Table I it is clear 
that single phonon emission cannot be invoked 
for an explanation. Although we are unable to 
resolve the puzzle, we should like to note here 
a significant difference between the valence band 
structures of Ge and Si which may be fruitful to 
explore further. This difference refers to the 
size of the spin-orbit splitting 4so of the valence 
bands at the central zone point. In Ge, on the 
one hand, Aco =0.29 > E boron) =0.01 ev. Hence 
the split-off valence band can play no direct role 
in the LFB of Ge. Further, the curvatures of 
the light- and heavy-hole bands remain substan- 
tially constant in the range 0 to 0.01 ev. In Si, 
on the other hand, a value of 0.05 ev has recently 
been reported® for Agog. Because of this, the 
curvature of the light-hole band changes abruptly 
at about 0.01 ev, and thereafter reverts toa 
value nearly the same as that for the heavy-hole 
band.” Thus, at energies comparable to E (boron) 
= Ago the effective-mass disparity between “light” 
and “heavy” holes has vanished, and therefore 
the interband scattering is changed in this energy 
region. Further, the “hot” holes may possibly 
be scattered into the split-off band before reach- 
ing energy E;. Further exploration is required 
to see whether these differences between Si and 
Ge can lead to energy-loss mechanisms which 
are much more effective in p-Si than in p-Ge.® 
It is of interest to note that the experimental 
evidence indicates that the light holes do not play 
an exceptional role in the LFB of p-Ge. Indeed 
the threshold rf power, at cyclotron resonance, 
for LFB in p-Ge is approximately the same for 
light and heavy holes.® 

Of the deep-lying, double acceptors in Ge, only 
Zn, according to currently available measure- 
ments (Table I), has a level with E;<E,.4. There- 
fore of this group of impurities (Zn, Cu, Ni, Fe, 
etc.), we expect only Zn to exhibit LFB in Ge. 
However, the lowest Cu level has an E; close 
enough to E opt that more accurate data on E; for 
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Cu would be desirable. The interpretation of dat, 
on breakdown produced by dc, or pulsed dec, fiej4, 
in Ge doped with the “deep-lying” acceptors mug 
proceed with great caution. First, care must bh. 
taken to establish that the observed breakdown j; 
not actually double-injection breakdown,” which 
yields a similar current-voltage characteristic, 
This latter breakdown is due to the simuitaneous 
injection into the sample, at some threshold yo}t- 
age, of electrons from the cathode and holes 
from the anode. Observations of this type of 
breakdown are usually sensitive to the contacts 
employed on the sample. Second, if stray light 
is shining continuously on the sample the observa 
breakdown may be due to impact ionization of 
shallower impurities which function as “deep 
traps” (at low enough temperatures), and hence 
are heavily populated by carriers, under illumi- 
nation. For the first of these reasons we believe 
that the interpretation of the breakdown data for 
Ni- and Cu-doped Ge presented by Sclar and 
Burstein’ is open to question. On the other hand, 
their breakdown data for Zn-doped Ge have been 
reproduced by one of us (MCS) and appear to be 
independent of the contacts employed. Thus it 
is almost certain that the observed breakdown 
of Zn-doped Ge samples is not due to double in- 
jection. Because of the magnitude of the thres- 
hold breakdown field (about 400 v/cm), the ob- 
served breakdown is also probably not due to 
impact ionization of acceptor impurities which 
are shallower than zinc. Therefore it appears 
that LFB has been observed in Zn-doped Ge. 

Further experiments are in order to clarify 
the entire picture of LFB. Perhaps the experi- 
ments on B-, Al-, and Ga-doped Si should be re- 
peated at microwave frequencies (to eliminate 
contact problems), and at higher rf powers, an 
with a controlled, optically excited, background 
of free holes. Even if low-field breakdown is no 
observed in B-doped Si, the breakdown must 
occur at some high field, and it would be instru 
tive to compare this threshold field with the cor- 
responding fields for Al-doped and Ga-doped Si. 
It would also be interesting to determine the thre 
hold breakdown fields for a series of Ge-Si alloy 
samples of varying composition but with identical 
boron (or Al or Zn or As) content. Such measurt) 
ments would provide a stringent test of the ideas 
we have discussed here, and may also add to ow 
present information" on the energy-band struc- 
ture and phonon spectra for these alloys. Brea 
down of Ge doped with “deep-lying” acceptors, 
Zn, Cu, Cd, Ni, etc., should be studied more 
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carefully than in the past, namely at microwave 
frequencies. Cu-doped Ge is particularly inter - 
esting because E (Cu) is so close to Ey (Table 1), 
We also note that the absence of LFB in Si makes 
possible the observation of space-charge-limited 
currents in Si at low temperatures, an experi- 
ment which directly measures the amount of com- 
pensation in a semiconductor.’” 

As the technology of the III-V compounds im- 
proves, we may expect to obtain additional use- 
ful data on LFB. This breakdown has already 
been observed in n-InP,” p-InSb,"* and n-InSb,”* 
on samples doped with very shallow impurities. 

The authors are indebted to many of their col- 
leagues for critical comments on this work. They 
would also like to thank Professor E. Burstein 
for several illuminating conversations on LFB. 





‘Electrical breakdown resulting from impact ioniza- 
tion across the band gap can be produced in any mate- 
rial at a sufficiently high field strength &,,.. In this 
note we are concerned with threshold breakdown fields 
$urp“SpG- For Ge and Si, &§ BG~ 10° v/cm and 
therefore the threshold breakdown fields under discus- 
sion here are ; pp <10* v/cm. For a list of refer- 
ences on high-field breakdown (across the band gap) of 
Ge and Si, see J. B. Gunn, Progress in Semiconduc- 
tors (John Wiley and Sons, New York, 1957), Vol. 2, 
p. 213. 

'LFB in Ge was first reported by E. Burstein et al., 
Phys. Rev. 91, 215 (1953); 93, 65 (1954). 

'N. Sclar and E. Burstein, J. Phys. Chem. Solids 
2, 1 (1957). 

‘§, H. Koenig and G. R. Gunther-Mohr, J. Phys. 
Chem. Solids 2, 268 (1957). 

‘Breakdown was sought, and not found, in many dif- 
ferent samples of n- and p-Si at 4°K with applied dc 
fields of up to 1000 v/cm, by one of the authors (MCS). 
R. D. Larrabee, at this Laboratory, was unable to 
produce breakdown in boron-doped Si at 4°K applying 





a pulsed dc field up to 60 kv/cm. The University of 
California (Berkeley) and Lincoln Laboratory groups 
studying cyclotron resonance in solids were unable to 
ionize shallow impurities in Si with the rf fields (in a 
microwave resonant cavity) available from a low-power 
klystron: Dexter, Lax, Kip, and Dresselhaus, Phys. 
Rev. 96, 222 (1954). 

*L. Huldt and T. Staflin, Phys. Rev. Lett. 1, 313 
(1958). The theoretical estimate for Aco is also close 
to this value: W. Kohn, Advances in Solid State Physics 
(Academic Press, New York, 1957), Vol. 5, p. 298. 

'E. O. Kane, J. Phys. Chem, Solids 1, 83 (1956). 

*another possibility which should not be ruled out is 
that the inelastic excitations of the boron acceptor 
states are considerably more effective in Si than in Ge 
for cooling “hot” holes. Note that the acceptor level 
structure in Si is different from that in Ge, i.e., not 
obtainable simply through a scaling factor, for two 
reasons. First, because of the smallness of Aco in 
Si, there is greater mixing-in of split-off band Bloch 
functions into the acceptor state wave function. Second, 
because of the greater depth of the acceptor ground 
states in Si, there are greater “core corrections’ to 
the effective-mass calculation. See W. Kohn, reference 
6, pp. 289 and 303. 

°R. C. Williams (private communication). This re- 
sult is consistent with strong scattering of the light 
holes into the heavy-hole band. This is also suggested 
by the almost identical fractional line-widths of the 
light- and heavy-hole cyclotron resonances [Dressel-~ 
haus, Kip, and Kittel, Phys. Rev. 98, 368 (1955)]. 

See also Goldberg, Adams, and Davis, Phys. Rev. 
105, 865 (1957). 

uble-injection breakdown has been reported for 
Fe-doped Ge by W. W. Tyler, Phys. Rev. 96, 226 
(1954), and for Cu-doped Ge by P. J. van Heerden, 
Phys. Rev. 108, 230 (1957). 
‘Braunstein, Moore, and Herman, Phys. Rev. 109, 
695 (1958). Herman, Glicksman, and Parmenter, 
Progress in Semiconductors (Haywood and Company, 








London, 1957), Vol. 2, p. 3. 

‘tM. A. Lampert, Phys. Rev. 103, 1648 (1956). 
13M. C. Steele, J. Phys. Chem. Solids 9, 93 (1959). 
‘R. Bray, Phys. Rev. 100, 1262(A) (1955). 

RR. J. Sladek, J. Phys. Chem. Solids 5, 157 (1958). 
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MEASUREMENT OF DIELECTRIC SURFACE CHARGE DENSITY 


Henry Seiwatz 
Physics Division, Armour Research Foundation, Chicago, Illinois 
(Received February 24, 1959) 


In a recent Letter, Shao’ reported a method for 


determining the quantity of charge transported 
across an electrode-dielectric interface during 
the charging of a capacitor. His technique is to 
measure the surface charge density which re- 
mains on the dielectric after the removal of one 
electrode of the charged capacitor. He assumes 
that the dielectric surface charge density of the 
intact capacitor is unaltered by this procedure. 
The purpose of this Letter is to point out that, 
except at low voltages, Townsend discharge oc - 
curs in the air gap as the electrode is removed, 
thus producing a considerable change in the sur- 
face charge density of the dielectric. As shown 
below, the dielectric surface charge densities 
calculated on the basis of Townsend discharge 
alone are of approximately the same magnitude 
as the experimental values reported by Shao. It 
is concluded that, in the proposed method of 
measurement, Townsend discharge masks the 
charge transfer which occurs across the elec- 
trode -dielectric interface of an intact capacitor. 


The minimum dielectric surface charge density 


which is produced by Townsend discharge alone 
can be calculated readily for a plane parallel 
capacitor. Let a capacitor charged to voltage V 
have a dielectric of thickness d and dielectric 
constant € as well as an air gap of length x. If 
the difference of potential across the air gap is 
limited by the sparking potential V, of air at 
atmospheric pressure, a surface density of 
charge o is formed on the dielectric. Gauss’ law 
gives 


V = 0d/e€,+(V./x)(x+d/e), (1) 


where ¢, is the permittivity of free space. When 
x is varied from zero to infinity, the surface 
charge density remaining on the dielectric is 


Oy, = (€€o/d)(V - V4), (2) 


where V; is the minimum of (V,/x)(x+d/e) with 
respect to x. 
This calculation can be applied to the geometry 


used by Shao. His measurements on mica (refer- 


ence 1, Fig. 2) are represented by graphs of the 
total charge density 2, of the capacitor and the 
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dielectric surface charge density 2 as functions 
of the voltage V over the range 0-15 kv. The 
effective capacitance C per unit area of the die- 
lectric is defined as = t/ V. For voltages in the 
range 0-7 kv, C is a constant equal to 140 
uuf/cm*. Above 7 kv, C increases with in- 
creasing voltage. The increase in C at higher 
voltages can probably be attributed, at least in 
part, to charging of the dielectric surface out- 
side the perimeter of the electrode by Townsend 
discharge. The constant value of C at low volt- 
ages is used to obtain an approximate value of 
d/e =6.3 microns for the mica dielectric. With 
this value of d/e and empirical values for the 
sparking potential’ in air, V; is found to be 0.54 
kv. Equation (2) then gives 


Oy, =1.4x10"% (Vv - 540), (3) 


for the voltage range in which C is constant, 
where V is in volts and o,,, is in coul/cm’. 

The experimental 2, vs V data can be fitted by 
an equation of the form 


Zq=C(V-V,), () 


in the voltage range 3-7 kv with V, =0.8 kv and 
C=140 ppf/cm*. The agreement between o,, 
calculated from Eq. (2) and the empirical £, is 
reasonably good in this voltage range. A meanin- 
ful value of o,,, cannot be calculated from Eq. (2) 
for voltages above 7 kv, because the mechanism 
responsible for the increase in C in this voltage 
range cannot be specified with assurance. A 
quantitative comparison of o,, with D4 at higher 
voltages is therefore not possible, but it is clear 
that they are of the same order of magnitude. 
Equation (2) is not in agreement with the Dy vs 
data at 1 kv and 2 kv, suggesting that electrode- 
dielectric charge transfer by mechanisms other 
than Townsend discharge is of considerable im- 
portance at very low voltages. 





‘Pp, Shao, Phys. Rev. Lett. 2, 41 (1959). 

23. M. Meek and J. D. Craggs, Electrical Breakdow 
of Gases (Oxford University Press, London, 1953), 
p. 84. 
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ANALYSIS OF DEUTERON STRIPPING EXPERIMENTS 


R. D. Amado 
University of Pennsylvania, Philadelphia, Pennsylvania 
(Received April 7, 1959) 


Deuteron stripping experiments analyzed accord- 
ing to the theory of Butler! are a nearly unique 
source of information on the orbital angular mo- 
mentum and single-particle widths of nuclear 
pound states. Although a vast amount of such 
information has been obtained in recent years, a 
number of problems in the understanding of the 
Butler theory remain. In particular the orbital 
angular moment of the captured particle can nor- 
mally be extracted unambiguously from the ex- 
perimental data, but the reduced width cannot. 
Furthermore the reasons for the success of the 
Butler theory or for its failures, at large angles 
for example, are not clearly understood. 

We wish to point out that recent work by Chew 
and Low in a different context cast light on these 





questions.* They show that in reactions in which 
there is a contribution from the exchange of a 
single particle, there can appear isolated poles 
in the renormalized Born approximation to the 
cross section, usually for unphysical scattering 
angle, and that the residue at these poles can be 
simply related to quantities of physical interest.® 
Stripping is such a reaction and the Butler theory 
is the renormalized Born approximation. To see 
this let us write the Butler theory for the process 
A(d,p)B. The Feynman graph corresponding to 
the lowest order Born approximation is shown in 
Fig. 1. The cross section in the center of mass 
corresponding to this diagram with the vertices 
and propagators treated exactly, that is renor- 
malized, is’ 





3 2 

do Os, +0) ame, 2M , a 1 6 
in - 2\2 2 
dQ (27, +1) Ri Mp 1 ar, (q’ +k ) (27 +1) 


x) YG, _ (GR) - (1+1)9, , 1 (GR) - tk 7,(QR) 


where R is the nuclear reaction radius, 7; the 
triplet effective range, and a is related to the 
deuteron binding energy, By: a*h?/m=Bg, where 
mis the nucleon mass. k,i is the wave number 

of the captured neutron, captured into a state 

with orbital angular momentum |. 6 is related 
tothe reduced width y by 7” =3h76?/2mR. My, and 
Mp are the masses of the target and daughter 


’ nuclei. The term in brackets contains the usual 


spherical Bessel functions. All wave numbers 
are appropriate to center-of-mass motion. 

The terms in Eq. (1) are easily interpreted in 
terms of Fig. 1. The factor a/(1- ar,) is related 
to the normalization of the asymptotic state of 
the deuteron and hence is the probability for dis- 
association corresponding to vertex 1.* The en- 
ergy denominator or propagator for the neutron 
intermediate state is 1/(¢* +k,,”). At vertex 2a 
neutron with orbital angular momentum / must 
be captured at radius R. The probability that 


F . 2 
lh (ik R) ~(1+1)h) (ck R) 


I-1 
- , (1) 
h (ik R) 








this partial wave and its derivative are in the 
neutron plane wave with momentum q is given 
by the term in brackets, and the probability that 
this wave will be captured is just the reduced 
width, ©. The other factors are kinematical. 
Equation (1) has a second order pole at the un- 
physical momentum transfer ¢ =-k,,”. Using 
energy conservation, we see that this corre- 
sponds to a scattering angle with cos@>1, and 
hence to an angle beyond zero, but for large in- 
cident deuteron energy, not very far beyond. At 
the pole all other contributions to the stripping 
process, for example compound nucleus, etc., 
are finite. Thus at the pole the Butler theory is 
exact, since it is infinitely larger than all other 
contributions, and the residue at the pole is ex- 
actly given by Eq. (1). If one measures the cross 
section, the only unknown in Eq. (1) is the re- 
duced width. This may be found by finding the 
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residue at the pole, that is by dividing the exper- 
imental cross section by the theoretical one, ex- 





do Ry (2J,+1) /M 


A B 





(27, + 1)6* = 


aQqk 3R 
p 


where now do/d&. is the experimental cross sec- 
tion. Figure 2 presents the results of such an 


analysis for the reaction C’*(d,p)C**.5 We see that 
the data allow an unambiguous extrapolation to the 


pole. 
In our analysis we exploit the dependence of 


the Born term on momentum transfer through the 
energy denominator 1/(¢ +k,,”). Usually it is the 
Bessel function term that is of interest since it 
has a strong dependence on /. For extracting the 
reduced width, this term can be a hindrance 
rather than a help once / is determined since it 


can vanish, or be small at small angles, and thus 


can give trouble when we divide through by it and 
try to extrapolate. Figure 2 represents a par- 
ticularly favorable case, but it is by no means 
unique. A further complication can be caused at 
very small angles by the Coulomb scattering. 
This is partly avoided by using high incident 
deuteron energies and low-Z targets, but even 
for higher Z it should be possible to make the 





FIG. 1. Lowest order Feynman graph for the reac- 
tion A@, p)B with incoming deuteron wave number kj, 
outgoing proton Kp, and momentum transfer q. 
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2M , 





cept for the reduced width, and then extrapolating 
to the pole. Doing this, we obtain 


l-ar, (v+k 7? 
t n 
{ P G=-k 2 


n 





(21+1) ; (2) 





extrapolation if one uses stripping theory with 
Coulomb corrections, so long as the Coulomb 
field can be treated adiabatically,® since then 
the particles may be thought to “ride up” the 
Coulomb barrier and “slide down” the other side 
without polarization of the internal deuteron mo- 
tion. So long as this is so the Born approxima- 
tion, with Coulomb wave functions,’ will still be 
very large at q =-k,,” and the extrapolation may 
be used. 

From our point of view the validity of the But- 
ler theory rests in the near vanishing of the en- 
ergy denominator at small angles. Thus the 
Butler theory, or Born approximation in general, 
works not only when interactions are weak, but 
also when energy denominators are small. In 
fact if the energy denominator is small, the Born 
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FIG. 2. The function of Eq. (2) plotted against the 
cosine of the center-of-mass scattering angle for the 
reaction C!(d, p)C"8, C'® being left in the 3. 08-Mev 
excited state,® with /=0 and R=4.2f. Extrapolation 
to the pole at cos@= 1.15 gives (2J9+ 1)67=12.5% 
The experimental differential cross section is also 
shown plotted on the displaced right-hand scale. 
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approximation will get better the larger the in- 
teraction matrix elements. At larger angles the 
Butler theory fails because it is no longer “riding 
the shoulder” of the pole. 

It is clear that methods such as these may be 
used in the analysis of a large number of direct 


interactions in which a single particle is exchanged. 





‘supported in part by the National Science Foundation. 
1a comprehensive list of references may be found in 
s. T. Butler, Nuclear Stripping Reactions (John Wiley 





and Sons, New York, 1957). 

2G. F. Chew and F. E. Low, Phys. Rev. (to be pub- 
lished). 

5The analytic properties of the scattering amplitudes 
have not been investigated rigorously. We shall assume 
that the extrapolation to the pole is permitted. 

‘J. M. Blatt and V. F. Weisskopf, Theoretical estene 
Physics (John Wiley and Sons, New York, 1952), p. 611. 

‘Freemantle, Gibson, and Rotblat, Phil. Mag. 45, 
1200 (1954). 

6s. T. Butler, reference 1, p. 62. 

"See W. Tobocman and M. H. Kalos, Phys. Rev. 97, 
132 (1955). 








EFFECTS OF ATOMIC ELECTRONS ON p-p AND n-p SCATTERING. I* 


G. Breit 


Yale University, New Haven, Connecticut 
(Received April 7, 1959) 


It has been pointed out’ that Coulomb excitation 
of atomic electrons by the incident proton and the 
excitation of atomic electrons caused by the accel- 
eration of the target proton in a proton or neu- 
tron collision may be expected to take place with 
appreciable probability. It was also brought out 
that in classical mechanics these effects do not 
matter because of the smallness of the energy 
changes and that in a quantum treatment the 
classical phenomenon of the disappearance of an 
electron with original energy and its reappearance 
at nearly the same angle with slightly changed 
energy is replaced by the effect of cross product 
terms on coherent scattering with the incident 
wave. For special channels these effects did not 
compensate’ the incoherent scattering. It ap- 
peared possible, therefore, that corrections for 
atomic electron effects to low-energy p-p scat- 
tering comparable with those for vacuum polari- 
zation come under consideration. On the other 
hand, the existence of compensations at high en- 
ergies due to the smallness of the proton wave- 
length in comparison with atomic dimensions was 
realized. Calculations by de Wit, Fischer, and 
lickendraht? have since then shown in the cases 
of monopole, dipole, and quadrupole excitation 
exact compensation of coherent scattering inter - 
ference effects by inelastic scattering provided 
contributions of all channels for different orbital 
iugular momenta L of the proton wave are com- 
bined and the following approximations are made: 
4) the asymptotic phase kr - $11 - nln(2kr) 
+argl'(L +1+4%n) of the regular Coulomb function 
is used and a corresponding replacement of the 


true by the asymptotic phase is made for the 
s-wave scattering anomaly, (b) in the evaluation 
of radial integrals all high-frequency parts in- 
volving 2kr in the phase are neglected. This 
compensation was interpreted as an indication of 
the smallness of effects of electron Coulomb ex- 
citation on observed scattering. It proved diffi- 
cult, however, to extend the method of summation 
over L to the evaluation of the residual effect, 
improvements on the asymptotic phase leading to 
difficult sums. On the other hand, a general 
reason for the compensation becomes apparent 
as previously expected’ for high-energy scatter - 
ing, the nucleon wavelength being short in com- 
parison with atomic dimensions even at + Mev. 
The present note is based on calculations employ- 
ing this viewpoint. Its formal aspects had been 
previously carried out for Coulomb excitation of 
nuclei® employing classical action functions S 
and S,, for the ground and nth excited states. 
These S give solutions of the wave equation by an 
extension of the exp(iS/#) representation of the 
wave function for the case of one S. Classical 
mechanics for relative motion with atomic quan- 
tum excitation is obtained by keeping the first 
two terms in the expansion of each S in powers 
of %. This approximation, for which the com- 
pensation effect is practically exact, is justifiable 
provided the dimensions of the region occupied 
by the perturbing energy are large in compari- 
son with the nucleon wavelength, a condition 
satisfied for the hydrogen molecule as well as 
the atom. Within wide limits the equivalence to 
classical mechanics holds independently of the 
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degree of adiabaticity of the excitation process. 
For the atomic case and Coulomb scattering 
combined with nuclear phase shifts, the arrange- 
ment of the calculation in terms of S, and S,, has 
been related to the expression for the Coulomb 
modification of the plane wave ~° so as to make 
estimates of quantum corrections to classical 
effects possible. The asymptotic form of ~° for 
kr >1, including correction terms to the plane- 
wave and outgoing-wave parts 0(1/kr) and cor- 
responding approximations to the three-dimen- 
sional Green’s function, give for any multipole 
exact compensation without the 0(1/kr) correc- 
tions. With their inclusion there appear effects 
involving an extra power of 1/(kay), ay =0.529 
x10 cm, which at 1 Mev are of the general 
order of 10~* of the Coulomb scattering and are 





negligible. A more detailed report is in prepara- 
tion, but the appearance of an extra factor 1/(kay) 
removes doubts regarding low-energy nucleon- 
nucleon scattering as a method for dealing with 
nuclear forces. 





"This research was supported by the U. S. Air Force 
under a contract monitored by the Air Force Office of 
Scientific Research of the Air Research and Develop- 
ment Command. 

1G. Breit, Phys. Rev. Lett. 1, 200 (1958). 

*de Wit, Fischer, and Zickendraht, Proc. Natl. Acai, 
Sci. U. S. (to be published). The writer is indebted 
to de Wit, Fischer, and Zickendraht for discussions 
regarding their calculations. 

3G. Breit and R. L. Gluckstern, Encyclopedia of 
Physics (Springer-Verlag, Berlin-Gottingen-Heidelberg, 
to be published), 41;, p. 496, Sec. 2. 








EMISSION OF ALPHA PARTICLES FROM A COMPOUND SYSTEM 
OF HIGH ANGULAR MOMENTUM* 


W. J. Knox, A. R. Quinton, and C. E. Anderson 
Yale University, New Haven, Connecticut 
(Received April 2, 1959) 


The energy spectra of alpha particles emitted 
following bombardment of Ni with high-energy 
oxygen ions have been measured. Emitted parti- 
cles pass through a proportional counter and into 
a CsI crystal. The energy of a particle is meas- 
ured in the CsI crystal while its dE /dx is meas- 
ured in the proportional counter. The relation 
between E and dE /dx identifies the particle. 
Spectra are recorded on a 20-channel analyzer 
gated by coincidences between E and dE /dx pulses 
of selected amplitudes or are observed and photo- 
graphed on an oscilloscope as a display of dE /dx 
vs E for individual particles. 

The spectra of alpha particles emitted at var- 
ious laboratory angles from the bombardment of 
a Ni target with 162-Mev O ions are shown in 
Fig. 1. Observations were made at laboratory 
angles of 17°, 31°, 46°, 60°, 90°, 110°, 134°, and 
158°. For clarity, only three of the spectra are 
shown in Fig. 1. The peaks of all of the spectra 
are shifted toward progressively lower energies 
as the angle of observation is increased. The 
total number of particles in any spectrum de- 
creases rapidly with increasing angle up to about 
90 degrees and then remains approximately con- 
stant at backward laboratory angles. The shapes 
of the spectra indicate that the total numbers of 


402 


alpha particles would actually increase at back- 
ward angles were it not for the effect of the low- 
energy cutoff of the detecting system. 

In Fig. 2 are shown the spectra transformed 
into the center-of-mass system of a 162-Mev 
O ion incident on a Ni target nucleus. Again for 
clarity only a few of the spectra observed are 
shown. Laboratory angles of observation are 
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FIG. 1. Energy spectra in the laboratory system of 
alpha particles emitted in the bombardment of Ni by 
162-Mev oxygen ions. Laboratory angles of observa- 
tion are indicated. 
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FIG. 2. Energy spectra in the center-of-mass sys- 
tem of emitted alpha particles. The ordinate is pro- 
portional to differential cross section per Mev per 
steradian. Laboratory angles of observation are in- 
dicated although each laboratory angle corresponds to 
various center-of-mass angles (see text). 


indicated. It should be noted, however, that these 
are continuous spectra and that particles of dif- 
ferent energies observed at a given angle in the 
laboratory come from different angles in the 
center-of-mass system. For example, particles 
in the laboratory spectrum at 90° correspond to 
center-of-mass angles of 119° at 8 Mev (center- 
of-mass energy), and 105° at 25 Mev. The spec- 
tra in the center-of-mass system all peak at 
approximately the Coulomb barrier energy (12 
Mev) to be expected for the emission of alpha 
particles from the system O plus Ni. The rela- 
tive differential cross sections are now at a 
minimum at 90° center-of-mass angle and strong- 
ly peaked in the forward and backward directions. 
The spectra appear to be predominantly of the 
nuclear evaporation type’ except for a pronounced 
forward peaking of the alphas at about twice the 
peak energy in the center-of-mass system. Pro- 
ton spectra from this reaction show similar 
features. 

Figure 3 shows the angular distributions in the 
center-of-mass system of alpha particles of var- 
lous energies. The values for the relative. cross 


sections have been taken from spectra plotted as 

in Fig. 2. Here it is seen more clearly that most 
of the alpha particles fall into distributions approx- 
imately symmetric about 90° with maxima in the 
forward and backward directions. There is some 
indication of forward peaking at small angles at 

all energies and the spectra are predominantly 
peaked in the forward direction at the highest 
energies. 

The angular momentum involved in the forma- 
tion of a compound system by a 162-Mev O nu- 
cleus incident upon a Ni nucleus has a mean ef- 
fective value of about 50 in units of f. Particles 
from such a system should be preferentially 
emitted in directions such that their orbital an- 
gular momenta are aligned with the initial angu- 
lar momentum of the system.?»* In the limit of 
high angular momentum and decreasing density 
of levels with increasing angular momentum, an 
angular distribution of particles proportional to 
(siné)~* is expected.* The distributions in Fig. 3 
are consistent with (siné)~’ in the backward di- 
rection for the angles observed and rise more 
rapidly than (siné)~* in the forward direction. 

We interpret the spectra as indicating that this 
reaction proceeds primarily by the formation of 
a compound system of high excitation from which 
alpha particles, among other particles, are 
emitted. The angular distributions of the alpha 
particles show the effects of the extremely high 
angular momenta with which these systems are 
formed. The asymmetric portions of the spectra 
presumably indicate the presence of some direct- 
interaction process or processes. It is perhaps 
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FIG. 3. Angular distributions of emitted alpha 
particles in the center-of-mass system. The ordinate 
is proportional to differential cross section per Mev 
per steradian. The dashed curve in the center figure 
is proportional to (sin@)~', the predicted distribution 
in the limit of high angular momentum. ? 
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significant that the most prominent feature of the 
forward -peaked alpha particles is their relatively 
large number at about 40 Mev in the laboratory 
system. Alpha particles of 40 Mev have the same 
velocity as that of the incident oxygen ions. 

The symmetric parts of the energy spectra in 
the center-of-mass system will be compared to 
nuclear evaporation calculations.’ Similar 
measurements on other reactions of this type 
are in progress. A detailed description of the 


——_ 


apparatus used will be published elsewhere. 





This work has been supported by the U. S. Atomic 
Energy Commission. 

'Dostrovsky, Rabinowitz, and Bivins, Phys. Rey. 
111, 1659 (1958). 

*L. Wolfenstein, Phys. Rev. 82, 690 (1951). 

8T, Ericson and V. Strutinski, Nuclear Phys. 8, 184 
(1958). 





X-RAY YIELDS IN ,.-MESONIC ATOMIC TRANSITIONS 


Jeremy Bernstein and Ta-You Wu 
The Institute for Advanced Study, Princeton, New Jersey 
(Received April 7, 1959) 


In a series of very careful experiments Stearns 
and Stearns’ have shown that the total K x-ray 
yield per stopped muon for muons stopping in 
various light materials is smaller by a consider- 
able amount than would be predicted on the 
grounds of any simple theory of the effect. In 
the actual experimental arrangement used by 
Stearns and Stearns the incident y-meson flux is 
known and the outgoing x-ray energy is meas- 
ured. The data are interpreted as follows: On 
the basis of a plausible picture of the capture 
process the » mesons, in the main, pass through 
the 2P u-mesonic atomic state. In this state the 
fate of the 1 meson is decided in a competition 
between the radiative 2P—1S transition and a 
2P—1S transition in which the energy released 
is taken up by an Auger electron ejected into the 
continuum. Other nonradiative processes such 
as direct muon capture by the nucleus contribute 
negligible branching ratios. Since the Auger 
probability is roughly independent of Z and since 
the radiative transition goes like Z‘*, we might 
expect, on this picture, that the x-ray yield Y 
would go like 


Y =const Z*/(C, +Z*const), 


where C, gives the Auger probability. Indeed 
Stearns and Stearns find that empirically the 
yield as a function of Z has this form, but that 
the constant C, determined from experiment 
would require an Auger probability some 300 
times larger than that computed by standard per- 
turbation-theory methods.” In an attempt to ac- 
count for this large discrepancy, Day and Morri- 
son® have invoked collision mechanisms involving 
the .-mesonic atom and its neighbors in the 
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dense material. Since this explanation in terms 
of collisional de-excitation has gained consider- 
able acceptance,‘ we felt that it would be worth- 
while to indicate why, in our belief, it cannot be 
correct. 

According to their argument the yp. -mesonic 
atom may undergo a collision with a neighboring 
atom whose electrons take up the 2P—1S de- 
excitation energy. Day and Morrison estimate 
that the cross section for this, o(2P~1S)=74,", 
where a, is the electron Bohr radius. The fact 
that the collision cross section emerging from 
their calculation depends on a, and not on a w 
the muon Bohr radius, is what gives them the 
factor of 300 needed to explain the experimental 
results of Stearns and Stearns. In this discus- 
sion the usual perturbation theory calculation of 
the Auger effect is not questioned. Now the 
collision de-excitation process of Day and Morri- 
son may be looked at as an “external” Auger ef- 
fect, i.e., one taking place during the time of 
collision ¢. It can readily be seen that the cross 
section 7a,” implies a transition probability of 
the order 10'*/sec for this “external” Auger ef- 
fect; whereas the usual perturbation theory gives 
the value ~10"'/sec for the same Auger transi- 
tion, 2P~1S, in a free atom. It appears to us 
that it is a priori unlikely that the Auger process 
which is inefficient for taking up the 2P—1S non- 
radiative transition energy in an isolated atom 
could be made so much more efficient by colli- 
sion with another atom. More exactly speaking, 
we asked why should an electron in a neighboring 
atom be ejected by an Auger process ~ 300 times 
more efficiently than an electron in the given 
atom itself? In fact, if we estimate the cross 
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section of the de-excitation of the 1 -mesonic 
atom from the 2P to the 1S state by a colliding 








electron, the Born approximation gives 


o(2P - 1S) =2n Jo(@)dcosé, (1) 
where 


2°z° 





K {mY , 
—- (=) *u [at eKe 


K, and K being the initial and final momenta of 
the electron. Thus 
2 

o<ma,. 
Because of the factor of a,” which appears, this 
collision de-excitation cross section is much too 
small to play any role in explaining the discre- 
pancy. In fact, any collision mechanism® in- 
volving the transfer of the large energy of the 
2P-1S transition of the .-mesonic atom to an 
electron will have a very small cross section on 
account of the great disparity in the masses of 
the meson and the electron. Hence, in our opi- 
nion, the experiment of Stearns and Stearns re- 
mains unexplained. We do not know of a mecha- 
nism which will account for the data. However, 
it is our belief that the usual perturbation-theory 
calculation of the Auger effect cannot be trusted. 
The reason is that the Auger perturbation and 
the other potentials in the Hamiltonian are, es- 
pecially for light elements, of the same order 
of magnitude and any separation of the Hamilto- 
nian is arbitrary. One of us* has pointed out 
that widths of atomic levels subject to the Auger 
effect are also very poorly estimated by the per- 
turbation theory. 

In order to tie the discrepancy down more 
sharply, it would seem very important to do an 
experiment in which the Auger rate is measured 
directly. A possibility’ is to do a coincidence 
measurement between 3D—2P or 3S—2P photons 
produced by the cascading muon and the subse- 
quent Auger electrons. This would measure the 
combined “external” and “internal” Auger rates. 
A material like carbon would seem particularly 
suitable. In the experiments of Stearns and 
Stearns, the Auger electrons were not actually 





-2KK, cos0?* a? + (9/4)Z"]°’ 


(2) 





seen. 

Then, in summary, the discrepancy between 
the experiment of Stearns and Stearns and the 
theory of x-ray yields seems to us to be unre- 
solved. 

We would like to thank Professor C. N. Yang, 
Professor L. Madansky, and Professor V. Te- 
legdi for helpful discussions and Professor J. R. 
Oppenheimer for the hospitality and support of 
the Institute for Advanced Study. 
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San example of another mechanism is an “exchange” 
Auger effect in which the muon makes a transition to a 
neighboring atom. In evaluating this process, there 
will be a factor corresponding to the electron efficiency 
in taking up the muon energy loss and a factor cor- 
responding to the probability of muon transfer. The 
muon can only make a transfer when the two atoms are 
within a distance a, /Z from each other, and hence this 
exchange process will have a cross section of at most 
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ONE-MESON CONTRIBUTION TO THE DEUTERON QUADRUPOLE MOMENT* 
David Y. Wong 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received April 6, 1959) 
According to Mandelstam’s representation, a’ =isin 26,61 . sind, ,- (2) 


partial-wave amplitudes for the nucleon-nucleon 
scattering have singularities only on the real 
axis of the momentum-square variable.’ These 
singularities correspond to physical processes 
such as the deuteron bound state, one or more 
meson exchange, etc. The low-energy triplet s- 
wave amplitude is quite accurately given by the 
deuteron pole alone. More precisely, the posi- 
tion and residue of the deuteron pole determines 
the scattering length and the effective range and 
vice versa. If one believes that the small mix- 
ture of d-state in the deuteron is mainly due to 
the long-range part of the nuclear force, then 
the residue of the deuteron pole for the mixing 
amplitude (mixing of s and d waves in the J=1 
scattering state) should be largely determined by 
the one-meson exchange contribution. When both 
the deuteron pole and the one-meson branch cut 
are included in the calculation of the mixing am- 
plitude, the requirement that this amplitude 
vanishes at zero momentum then gives the re- 
sidue of the deuteron pole in terms of the meson- 


nucleon coupling constant. This residue is simply 


related to the ratio of the asymptotic s- and d- 
wave functions of the deuteron. The result is in 
good agreement with those of the phenomenologi- 
cal deuteron wave functions which have been ad- 
justed to the observed quadrupole moment as 
well as other low-energy data. A conventional 
direct estimation of the quadrupole moment in 
terms of the asymptotic s-d wave ratio gives a 
value within 40% of the experimental result. 

In the notation of Stapp et al.,* the s-wave am- 
plitude a, , and the mixing amplitude a’ are 
given in terms of the eigenphases 59,1» & 1, and 
the mixing parameter e¢, by 


16 
Qo, 1 = 2i(cos*«e,e 9, : sind, 1 


+ sin*e,e°°™ * sind, ,), 


76 76, 
a‘=isin2e,(e sind, ,-e **sind, ,), (1) 


and a corresponding expression holds for the d- 
wave amplitude @,,,- For small ¢, and 61, Eq. 
(1) can be replaced by 


i6 
Qo, i= 2ie %, s sind, 1? 


In the effective-range approximation, we have 





@, 1 id, 1 

——=-—e »* sind, C eee 

2igq 91 qcot by , - ig 
- : (3 
* W/a)+i rg ig’ 


where q is the center-of-mass momentum of 
either one of the nucleons, a is the scattering 
length, and ry is the effective range. Experimen- 
tal values for a and ry in meson units are*® 


a=- 3.82, r=1.20. (4) 
It is clear from Eq. (3) that (a, ,/2iq) has poles 


at 
1 1 2\4 -1.79 
es -[+(4+2) ges 6) 


As anticipated, the value g,” = 0.106 is almost ex- 
actly the deuteron binding energy times the nu- 
cleon mass m (q*/m is the kinetic energy in the 
center-of-mass system). The pole at -1.79 can 
be interpreted as an interaction that gives rise 
to the bound state. The residue of the deuteron 
pole at -q,” is 


I = -2g,/(1 - qq) = - 1.07. (6) 


In the effective-range theory, -I is just the 
square of the normalization constant of the deu- 
teron radial s-wave function. To summarize, 
any two of the four parameters {a, r, q,”, I} de- 
termine the remaining two inasmuch as the deu- 
teron pole dominates the low-energy behavior of 
the a, , amplitude. 

Since the deuteron is known to have a small 
mixture of d-state, the mixing amplitude (a'/iq) 
will also have a pole at -g,”. From the nonrela- 
tivistic scattering theory,‘ the residue is (-2Ip), 
where p is the ratio of the asymptotic s- and d- 
wave functions of the deuteron. When the one- 
meson exchange process is included, a branch 
point will occur at q*=-4. If the branch cut is 
taken from - 4 to -«, the discontinuity of (a'/iq) 
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across the cut can be shown to be® 


BL) ae Ca) ae el 





q*+ie 
2inf?m /1 3 —— 
= 5 (= + iz") for q* <-%. (7) 


On the positive real axis of the qg* complex plane, 
the phase of (a*/ig) from Eq. (2) is given by 5), 
which can be calculated from Eq. (3) using the 
effective-range formula. Together with the pre- 
scribed singularities on the negative real axis, 

the solution of (a'/ig) can be found by the gen- 
eral method of Omnes.® However, for the present 
problem, it is simple enough to write down a 
solution that possesses all the required analytic 
properties. One obtains 


s tt a= -20p 
(at/iq) = +z are ia) Tae qe Golla*+a0) 





f?m “4 1/q*+3/4q" 
a ava aaaal 





This should be a reasonable approximation at low 
energy since it reproduces exactly the deuteron 
pole, the one-meson cut, and the low-energy 
phase shift. The requirement that (a'/ig) vani- 
shes at zero momentum then gives 


p= 0.029, (9) 
for f?=0.08. This value of p is consistent with 


those of the phenomenological deuteron wave func - 
tions which range between 0.025 and 0.032.” ® 

A direct estimation of the quadrupole moment 
in terms of p can also be made in the conventional 
way. This gives a quadrupole moment of ~ 4 
x 10°?" cm® compared with the experimental value 
of 2.71077 cm’. 

A more detailed consideration of the one-meson 
contribution to nucleon-nucleon scattering will be 
published in a separate report with H. P. Noyes. 

The author is indebted to Professor G. F. Chew 
for discussions which led to the present calcula- 
tion. 





"This work was done under the auspices of the U. S. 
Atomic Energy Commission. 

'S, Mandelstam, Phys. Rev. 112, 1344 (1958). 

Stapp, Ypsilantis, and Metropolis, Phys. Rev. 105, 
302 (1957). 

3For example, see H. A. Bethe and P. Morrison, 
Elementary Nuclear Theory (John Wiley and Sons, Inc., 
New York, 1956). 

‘For example, see Goldberger, Nambu, and Oehme, 
Ann. Phys. 2, 226 (1957). 

Scziffra, MacGregor, Moravcsik, and Stapp, Univer- 
sity of California Radiation Laboratory Report UCRL- 
8510, 1958 (unpublished). 

*R. Omnes, Nuovo cimento 8, 316 (1958). 

"L. Hulthen and M. Sugawara, in Handbuch der Physik 
(Springer-Verlag, Berlin, 1957), Vol. 39, p. 92. 

®Previous work on relating the coupling constant with 
the s-d wave ratio can be found in Iwadare, Otsuki, 
Tamagaki, and Watari, Progr. Theor. Phys. (Kyoto) 
15, 86 (1956), and 16, 455 (1956). 
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ABSTRACTS 








In this section are printed the abstracts of Articles that 
have been forwarded to The American Institute of Physics 
for publication in THE PHYSICAL REVIEW. In quoting 
information obtained from this section before the appear- 
ance of the corresponding Article, reference should be 
made to “Physical Review (to be published)” rather than 
to this Journal. 


BEHAVIOR OF INTENSE RELATIVISTIC ELEC- 
TRON BEAMS. Richard H. Capps,* University 
of Washington, Seattle, Washington (Received 
June 11, 1958). 


A theoretical investigation is made of intense, 
relativistic, partially neutralized, self-focusing 
electron beams. Approximate equations are de- 
rived that express the time derivatives of the 
beam radius and mean electron energy as func - 
tions of these parameters, the densities of the 
electrons, ions, and neutral atoms, and the im- 
pressed electric and magnetic fields. The equa- 
tions are numerically integrated for various 
values of the parameters. It is shown that the 
radiation damping leads to a shrinkage of the 
beam radius to a value of order 0.05 mm ina 
time of order 10™-10-' second. If the longitudi- 
nal electric field is constant the beam approaches 
the “equilibrium state” discussed by Budker, in 
which the increase in electron transverse mo- 
mentum caused by scattering from the ions is 
compensated by radiative loss, and the energy 
lost by radiation is restored by the longitudinal 
field. The various instability problems asso- 
ciated with dense systems of charged particles 
are not discussed. 


* 
Present address: Cornell University, Ithaca, New 
York. 


ANALYSIS OF THE LINEAR COOPERATIVE 
PROBLEM AS A MARKOFF PROCESS. W. W. 
Mullins, Metallurgy Department, Westinghouse 
Research Laboratories, Pittsburgh, Pennsylvania 
(Received November 17, 1958). 


The equilibrium statistical structure of a linear 
chain of interacting elements, each capable of 
a states, is determined by utilizing the theory 
of the entropy of a Markoff process to express 
the free energy of the chain as a function of its 
transition probabilities. Minimization of the 
free energy then leads to equations for the transi- 
tion probabilities. The condition that these equa- 


408 


tions have a nontrivial solution permits the de- 
termination of the free energy from the largest 
eigenvalue of a certain matrix (A), already famij- 
iar in the matrix method of cooperative phenom- 
ena. The corresponding eigenvector (together 
with the elements A;;) is shown to determine the 
statistical structure of the chain. Previous dis- 
cussions of the cooperative chain, in terms of 
a Markoff process, have been based on the 
assumption that the transition probabilities are 
given by normalized Boltzmann factors. This 
assumption is shown to be erroneous for most of 
the cases to which it has been applied. 


SINGLE-PARTICLE ENERGIES IN A DEGENER- 
ATE ELECTRON GAS AT HIGH DENSITY. Ronali 
M. Rockmore, Brookhaven National Laboratory, 

Upton, New York (Received December 11, 1958), 


Expressions for the self-energy of an electron 
and the physical pair-propagator for an excited 
electron plus hole in a degenerate electron gas a 
high density are obtained through the formulation 
of an equivalent Hamiltonian for the problem. 
These expressions are identical to those obtained 
by Quinn and Ferrell. A derivation of the physi- 
cal pair-propagator from ground state considera- 
tions is also given. 


HEAT CAPACITY OF ALUMINUM BETWEEN 
0.1°K AND 4.0°K. Norman E. Phillips, Depart- 
ment of Chemistry and Radiation Laboratory, 
University of California, Berkeley, California 
(Received December 1, 1958). 


Measurements of the heat capacity of alumi- 
num have been made between 0.11 and 4.0°K in 
the normal state and between 0.17 and 4.0°K in 
the superconducting state. Within the experimen- 
tal error the normal state heat capacity, Cy, cal 
be represented by C,, =yT+8T* with y =1.35x10" 
joule/mole deg? and a value of 8 corresponding 
to a Debye temperature of 427.7° in agreement 
with calculations based on elastic constants. For 
reduced temperatures between 0.5 and 0.25 the 
electronic heat capacity in the superconducting 
state, Cys, is approximated by C,,/yT,=7.1 
x exp(-1.34 T/T) in which T, is the transition 














temperature, 1.163°K. At reduced temperatures 
less than about 0.25, C,, is greater than an extra: 
polation of the exponential, the difference amoutt- 
ing to a factor of 4 at the lowest temperature. 

The departure of C,, from an exponential temper 
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ature dependence, which is believed to be outside 
the experimental error, is not consistent with the 
existence of a constant energy gap at low reduced 
temperatures. The calculated critical field is 
103.0 gauss at 0°K and shows a maximum nega- 
tive deviation of 4% from the parabolic law. The 
results are compared with other-measurements 
and with theory. 


CHARGE DISTRIBUTION AND RANGE EFFECTS 
PRODUCED BY 3-Mev ELECTRONS IN PLEXI- 
GLAS AND ALUMINUM. Bernhard Gross” and 
Kenneth A. Wright, High Voltage Research La- 
poratory, Department of Electrical Engineering, 
Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts (Received December 1, 
1958). 

Charge distribution curves of nearly parallel 
beams of monoenergetic electrons of 3 Mev have 
been measured in Plexiglas and aluminum. The 
charge distribution in depth reaches a maximum 
atapproximately 2/3 of maximum range. The 
maximum in Plexiglas is sharper than in alumi- 
mum. The results confirm the existence of a 
space-charge layer in dielectrics after electron 
bombardment. The measured charge distribution 
reflects the differential range distribution of 
monoenergetic electrons. 


‘on leave of absence from National Institute of Tech- 
nology, Rio de Janeiro, Brazil. 


ALTERNATIVE SUPERCONDUCTING GROUND 
STATE. J. C. Fisher, General Electric Research 
Laboratory, Schenectady, New York (Received 
December 15, 1959). 


A comparison is made between (a) the trial 
ground state composed of a linear combination of 
normal state configurations in which individual 
electron states of opposite momentum and equal 
spin projection (k+, -k#) and (k}, -k +), abbre- 
vated (ko, -ko), are both either occupied or un- 
occupied, and (b) the trial ground state with 
+, -k}) pairing proposed by Bardeen, Cooper 
and Schrieffer. The two trial states lead to iden- 
tical sets of one-electron matrix elements be- 
tween corresponding excited states, and to iden- 
tical thermodynamics as long as an average matrix 
element -V = (Ver ) ay is treated as a disposable 
parameter. They differ in that spin paramag- 
netism is predicted for (ko, -ko) pairing. De- 
tailed knowledge of the electron-phonon interac- 


tion in superconductors seems at present inade- 
quate to rule out either possibility for all super - 
conductors. 


IDENTIFICATION OF THE MAGNETOSTATIC 
MODES OF FERRIMAGNETIC RESONANT 
SPHERES. P. Fletcher, I. H. Solt, Jr., and 

R. Bell, Hughes Aircraft Company, Culver City, 
California (Received December 15, 1958). 


A detailed comparison between the magneto- 
static theory and experimental observation is 
given. Both the resonant field and the intensity 
of the magnetostatic modes are compared. The 
disagreement between observation of resonant 
fields and the static theory is about 50 gauss for 
an yttrium iron garnet sphere 1.3 mm in diameter. 
Inclusion of the first-order propagation correc- 
tions reduces this disagreement to less than two 
gauss for 29 of the 31 modes that were compared. 
This increased accuracy allows both a more 
positive identification and a more accurate deter - 
mination of g-factors. A comparison of the static 
theory of intensities with observation has like- 
wise been made. The observed intensities are 
close to the predicted intensities for some lines 
and from 4 to 100 times greater for others. 


OPTICAL PROPERTIES OF TELLURIUM AND 
SELENIUM. Robert S. Caldwell* and H. Y. Fan, 
Purdue University, Lafayette, Indiana (Received 
December 15, 1958). 


The optical properties of tellurium and of 
trigonal and amorphous selenium have been in- 
vestigated at wavelengths extending from the in- 
trinsic absorption edge to about 152 microns using 
polarized radiation. The refractive indices of 
tellurium and trigonal selenium have been deter- 
mined from 4 to 14 microns and from 9 to 23 mi- 
crons, respectively. For amorphous selenium, 
the refractive index estimated from measured 
reflectivity shows no appreciable variation from 
30 to 152 microns. The temperature shift and 
the pressure shift of the intrinsic absorption edge 
in amorphous selenium are found to be -1.45x1075 
ev/*K and -2.0x10~° ev/atmos, respectively. Lat- 
tice absorption bands have been observed in tellu- 
rium and selenium; they are attributed to the ex- 
citation of combination modes. In tellurium, a 
strong absorption band has been observed at 11 
microns, which is present only for E |ic radia- 
tion. The band indicates that there are overlap- 
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ping branches in the valence band which are sep- 
arated by about 0.11 ev. The effective mass in 
the lower branch is estimated to be about four 
times smaller than that in the upper branch. It 

is possible that the structure of the valence band 
is responsible for the high-temperature reversal 
of the Hall effect in tellurium. The usual carrier 
absorption increasing smoothly with wavelength 
has been studied for tellurium using polarized 
radiations. The effective mass of holes in tellu- 
rium has been determined from reflectivity meas- 
urements: m, ~m\ =0.45m at 300°K; m , =0.30m 
and m | =0.45m at 100°K. 


“Now at Boeing Airplane Company, Seattle, Washing- 
ton. 


APPLICATION OF THE MORSE POTENTIAL 
FUNCTION TO CUBIC METALS. L. A. Girifalco 
and V. G. Weizer, Lewis Research Center, Na- 
tional Aeronautics and Space Administration, 
Cleveland, Ohio (Received December 4, 1958). 


The Morse parameters were calculated using 
experimental values for the energy of vaporiza- 
tion, the lattice constant, and the compressibility. 
The equation of state and the elastic constants 
which were computed using the Morse parameters, 
agreed with experiment for both face-centered and 
body-centered cubic metals. All stability condi- 
tions were also satisfied for both the face-centered 
and the body-centered metals. This shows that 
the Morse function can be applied validly to prob- 
lems involving any type oi deformation of the cu- 
bic metals. 


COHERENT ATOMIC SCATTERING FACTORS 
FOR THE LITHIUM HYDRIDE CRYSTAL FIELD. 
R. P. Hurst, University of Wisconsin Naval Re- 
search Laboratory, Madison, Wisconsin (Re- 
ceived December 3, 1958). 


This paper describes an effort made to eval- 
uate the coherent atomic scattering factors for 
the lithium and hydride ions as they exist in the 
LiH crystal field. The two-electron wave func- 
tions chosen to describe these ions are each in 
turn optimized in the field of positive and nega- 
tive point charges arranged as the corresponding 
ions are arranged in the crystal. Using the or- 
bital exponents obtained from the minimization 
process, atomic scattering factor calculations 
are made in the usual way. The factors for Li* 
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are unchanged. Interesting effects are noted for 
the H~ scattering factors. 

The energy of the ions in the field of point 
charges is related to the crystal cohesive energy, 
Specifically, neglecting the zero-point energy 
correction, the calculated open configuration re- 
sults were -236.3 and -227.4 kcal/mole when 
compared with computed and experimentally de- 
termined free H~ reference energies, respec- 
tively. The experimental cohesive energy is 
-218.5 kcal/mole. 

Finally, it was noted that the H~ ion contracts 
in the crystal field. A possible explanation for 
this effect is given. 


TRANSPORT PROPERTIES OF NONDEGENER- 
ATE n-TYPE SEMICONDUCTORS CONSIDER- 
ING ELEC TRON-ELECTRON SCATTERING. 
M. S. Sodha,* Physics Department, University 
of British Columbia, Vancouver, Canada, and 
Y. P. Varshni, Physics Department, Allahabad 
University, Allahabad, India (Received Decem- 
ber 4, 1958). 


Sodha and Eastman have shown that the time of 
relaxation of an electron due to scattering by 
ions and electrons is given by 7=Bd(x)/x, where 
x is proportional to its velocity and ¢(x) isa 
function tabulated by Spitzer and Harm. In this 
paper the authors have proposed two analytical 
functions for ¢(x) which make the integrals ap- 
pearing in transport theory readily integrable. 


* 
Now at Physics Research Department, Armour Re- 
search Foundation, Chicago, Illinois. 


SURFACE POTENTIAL, FIELD EFFECT MOBIL- 
ITY, AND SURFACE CONDUCTIVITY OF Zn0 
CRYSTALS. H. J. Krusemeyer, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Received 
December 4, 1958). 


Measurements of the electrical properties of 
ZnO crystals were made at 300°K, both in dry 
O, and in high vacuum. The dark conductivity 
was changed by illumination with ultraviolet. Sur- 
face potentials of crystals with diameters down to 
0.002 cm were measured by the Kelvin method 
with a sensitivity of 0.002 v. Comparison of the 
measured and calculated dependence of dark con- 
ductivity on dark surface potential showed that the 
latter could be changed from about 0.1 v below to 
0.5 v above the neutral point. Application of a 
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transverse electric field produces a fast change 

in conductivity in less than 50 usec and a slow 
change in which part of the fast change decays 

with a time constant ranging from minutes to hours 
depending on ambient and surface potential. The 
field effect mobility increases with increasing sur- 
face potential from a value which is sometimes 
smaller than one to a plateau value between 70 and 
145 cm? v' sec™*, and decreases again for the 
largest value of surface potential. Evidence is 
given that the low mobility values are caused by 
surface states. Combined measurements of sur- 
face potential, field effect mobility, and surface 
conductivity together with quantum efficiency 
measurements of the surface conductivity by Col- 
lins and Thomas yield the quantum efficiency of 
the hole-trapping process at the surface which is 
approximately 1 for a neutral surface. A quanti- 
tative treatment of the hole-trapping process is 

in good agreement with the experimental results 
and shows that the bulk diffusion length for holes 
is >1000 A and that the ratio of hole surface 
trapping velocity and diffusion constant equals 
1.7x10° cm“. 


MOBILITY OF HYDROGEN IONS (H*, _. H,*) 
IN HYDROGEN. Edward A. Mason and Joseph T. 
Vanderslice, Institute for Molecular Physics, 
University of Maryland, College Park, Maryland 
(Received November 19, 1958). 


Force laws for H*, H,*, and H,* in H, are cal- 
culated from theory and from results on the scat- 
tering of low-velocity ion beams in H, gas. These 
results are then used to calculate the mobilities 
of the hydrogen ions in H, gas as a function of 
temperature. The mobilities of H+ and H,* de- 
crease slightly with increasing temperature, but 
the mobility of H,* increases strongly. The agree- 
ment with experiment indicates that the unidenti- 
fied hydrogen ion whose mobility has been meas- 
ured is probably H,*, rather than the usually 
assumed H,*. 


CLASSICAL THEORY OF ELECTRONIC AND 

IONIC INELASTIC COLLISIONS. Michal Gryzinski, 
Warsaw University, Hoza, Warsaw, Poland, and 
Institute of Nuclear Research, Warsaw, Poland 
(Received November 7, 1958). 


A classical theory of inelastic atomic collisions 
is given. It is shown that inelastic scattering, 
ionization, excitation, and other interactions 


between charged particles and atoms are due to 
the Coulomb interaction with atomic electrons 
and depend in a first approximation on their bind- 
ing energy and momentum distribution. All cross 
sections can easily be calculated by means of 
differential cross sections o(AE) and o(AE, 6) de- 
rived in the binary encounter approximation. Nu- 
merical calculations have been made for several 
cases and are in very good agreement with the 
experimental results. 


COMBINED USE OF THE METHODS OF SUPER- 
POSITION OF CONFIGURATIONS AND CORRE - 
LATION FACTOR ON THE GROUND STATES OF 
HELIUM-LIKE IONS. Per-Olov Lowdin and 
Lajos Rédei, Quantum Chemistry Group, Upp- 
sala University, Uppsala, Sweden (Received De- 
cember 1, 1958). 


The accurate solution of the Schrédinger equa- 
tion by means of a combination of the method of 
correlation factor and the method using super - 
position of configurations is discussed. Fora 
many-electron system, the total wave function 
divided by the nodeless function g=g(7,,, 715, 
Yo5---) May be expanded in a series of Slater de- 
terminants built up from a complete one-electron 
set. For a two-electron system, this expansion 
becomes very simple and, by going over to prin- 
cipal orbitals, it may be “diagonalized” and 
brought to a particularly rapidly convergent form. 

The method is applied to the He-like ions (H’, 
He, Li*, and Be?*) by choosing g=1+a7,, and by 
using a basis consisting of only three s-orbitals. 
The energies obtained differ from the exact val- 
ues only by about 0.001-0.002 atomic unit. Even 
the simpler wave functions (u)*(1+a7,,) and (u, v) 
x(1+a7,,) are discussed, and the forms obtained 
by truncating the expansion of ¥ in principal or- 
bitals are studied in greater detail. 


GROSS STRUCTURE IN (d,p) REACTIONS NEAR 
Z=50. R. A. Peck, Jr., Brown University, Prov- 
idence, Rhode Island, and J. Lowe, University 

of Birmingham, Birmingham, England (Received 
June 18, 1958). 


Proton spectra from Sn™(d,p)Sn'*! and In™*(d, 
p)In*** at 19.7 Mev are presented, covering 6.5 
Mev of excitation in the product nucleus and lab~ 
oratory angles from 15° to 155° at 10° intervals. 
The reactions are closely similar in all charac- 
teristics. Broad energy groups are discussed 
relative to the correlation with neutron single- 
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particle states anticipated by a giant-resonance 
interpretation of their nature. Of four such groups 
identified, the spectral locations of three agree 
quantitatively with the single-particle level scheme 
deduced from binding-energy data. A group at 

0.8 Mev, while consistent with single-particle 
gross structure, may represent a distinct nuclear 
level enhanced for reasons which are not clear. 
Angular distributions for the several groups fail 
to show the strong minima characteristic of pure 
stripping forms, probably as a consequence of 

the spreading of single-particle levels resulting 
from finite nuclear distortion which also weakens 
the spectral gross structure. However the slight 
structure which the distributions do display is 
qualitatively consistent with stripping-peak posi- 
tions consequent on the single-particle interpre- 
tation of the spectra. 


SEARCH FOR THE ISOTOPE Ir’. Donald G. 
Gardner* and W. Wayne Meinke, Department of 
Chemistry, University of Michigan, Ann Arbor, 
Michigan (Received December 8, 1958). 


The isotope Ir’® has been reported to have a 
half-life of ~9 days and emit 8~ particles with a 
maximum energy of about 0.08 Mev. Using deu- 
teron-bombarded enriched isotopes of platinum, 
it is shown that the previous mass assignment 
was incorrect. It is suggested that the ~ 9-day 
activity found in deuteron bombardments of nat- 
ural platinum is due to Ir*®* and Ir’™ produced by 
the (d,n) reaction on osmium impurities. An ex- 
perimental upper limit of 5 hours for the Ir’® 
half-life can be set by these experiments. Rough 
cross sections for the (d,a) reaction on Pt!™ and 
Pt are given for several deuteron energies 
from 9.6 to 20.4 Mev. 


* Present address: Westinghouse Electric Corporation, 
Pittsburgh, Pennsylvania. 


INELASTIC SCATTERING OF 21.6-Mev DEU- 
TERONS. J. L. Yntema and B. Zeidman, Ar- 
gonne National Laboratory, Lemont, [linois 
(Received October 16, 1958). 


The inelastic scattering of 21.6-Mev deuterons 
by Ni, Cu, Zn, Rh, Ag, Sn, Gd, Ta, Pt, and Au 
has been measured. A Nal(T1) crystal sandwich 
was used for the detection of the scattered par- 
ticles. The calculated energy scale of the spec- 
trum was checked by observing the inelastic 
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spectrum from Mg. The angular distribution 
shows a pronounced peaking in the forward di- 
rection. At backward angles there is no increase 
in cross section with increasing excitation, 
showing that compound nucleus processes do not 
contribute to the cross section. There is no 
measurable difference for the inelastic scatter- 
ing cross section for 21.6-Mev deuterons scat- 
tered at 60° by Au up to 10 Mev excitation as 
compared to 22-Mev protons. However, the deu- 
teron inelastic cross sections for Ni, Cu, and Zn 
are about half of the corresponding 22-Mev proton 
cross sections at 60°. The gross structure found 


by Cohen in the inelastic spectrum is also observed 


in the present experiment. However, the gross 
structure in Ta and Gd appears to be consider- 
ably more pronounced in the deuteron spectrum 
than in the proton spectrum. 


ELASTIC SCATTERING OF 21.6-Mev DEUTER- 
ONS BY SEPARATED ISOTOPES OF NICKEL 
AND COPPER. J. L. Yntema, Argonne National 
Laboratory, Lemont, Dlinois (Received Novem- 
ber 24, 1958). 


The relative differential cross sections for the 
elastic scattering of 21.6-Mev deuterons by Ni™, 
Ni®, Cu®*, and Cu® have been measured. The 
diffraction pattern shows a shift with atomic 
number. The shift from Ni®*® to Ni® is about the 
same as the one between Cu™ and Cu®. The 
shift between Ni® and Cu® is about 1.5 times 
that of the others. 


DEUTERON-PICKUP REACTION IN AN OPTI- 
CAL-MODEL APPROXIMATION. Kenneth R. 
Greider, Lawrence Radiation Laboratory, Uni- 
versity of California, Berkeley, California (Re- 
ceived December 12, 1958). 


A theory of the (p,d) pickup reaction is de- 
scribed in which the nuclear interactions of the 
incoming and outgoing particles are considered. 
Two different formal expressions that give the 
transition amplitude are derived, and the wave 
functions in this amplitude are approximated by 
an optical-model procedure in which it is as- 
sumed that the initial- and final-state particles 
scatter elastically in the nucleus. Several 
closed forms for these optical-model wave func- 
tions are derived on the basis of a WKB approxi- 
mation for a complex square-well scattering 
potential. The use of these wave functions, along 
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with an approximation that gives the form of the 
transition amplitude in terms of Gaussian func - 
tions, allows a closed-form solution for the 
differential cross section. 

It is found that the elastic -scattering processes 
are not negligible, since they affect considerably 
the magnitude and the shape of the differential 
cross section. By comparing the theory with re- 
cent pickup experiments on C™ at 95 and 145 
Mev, one obtains a nuclear-momentum distribu- 
tion that, unlike the Born approximation analysis, 
is in good agreement with the results of other 
determinations of momentum distributions. It is 
found that a neutron number of from 4 to 6 neu- 
trons and a momentum distribution of exp(-E/14) 
are required to fit the data. 


LIFETIME OF THE FIRST EXCITED STATE OF 
K®. F. J. Lynch and R. E. Holland, Argonne 
National Laboratory, Lemont, Illinois (Received 
December 8, 1958). 


The first excited state of K® at 29.4 kev has 
been produced by using the reaction A‘*°(p,n)K®. 
Its mean life, measured by a pulsed-beam tech - 
nique, is 5.6+0.5 mysec. 


RADIATIVE CAPTURE OF PROTONS IN Be’. 
J.J. Singh, Department of Physics, University 
of Kansas, Lawrence, Kansas (Received Decem- 
ber 15, 1958). 


The energy levels of B’° have been the subject 
of considerable experimental and theoretical 
investigations. In the present experiment, we 
have made some measurements on the capture 
gamma-rays in Be® at proton bombarding ener- 
gies of 662 kev and 998 kev. Branching ratios 
for various levels are given. 


NUCLEAR SPINS OF I’** AND ['*, J. E. Sherwood 
and S. J. Ovenshine, Oak Ridge National Labor- 
atory, Oak Ridge, Tennessee (Received Decem- 
ber 15, 1958). 


The nuclear spins of I'** and I’ have been 
measured by means of an atomic beam magnetic 
resonance apparatus. A value of 1 was obtained 
for I'*, as predicted by Benczer et al. However, 
a value of 5 was obtained for I’, in disagree- 
ment with the expectations of Caird and Mitchell. 
The nature of the discrepancy is indicated. 





The apparatus and technique are described, 
with special attention to the source and detector 
arrangements. 


LIFETIME OF THE FIRST EXCITED STATE OF 
Li’. C. P. Swann, V. K. Rasmussen, and F. R. 
Metzger, Bartol Research Foundation of the 
Franklin Institute, Swarthmore, Pennsylvania 
(Received December 15, 1958). 


The Li’(p, p’ )Li’* reaction was used as the 
source of Doppler-broadened 478-kev gamma rays 
to determine the mean life of the first excited 
state of Li’ by the nuclear resonance fluores- 
cence technique. Both scattering and self-ab- 
sorption experiments using both LiF and Li as 
the scattering and absorbing materials were 
performed. 

The resonance scattering gives a mean life of 
(1.09 +0.07) x10~* sec after 2% correction for 
slowing down of the Li”” nuclei has been made. 
This correction follows from the measured stop- 
ping cross sections for lithium ions. An addi- 
tional correction, for the effect of collisions in 
which the Li” nuclei lose a large fraction of 
their energy, is found difficult to estimate, al- 
though it is certainly not more than a few per- 
cent. For this reason, the scattering measure- 
ment is regarded as giving a lower limit for the 
mean life. 

The self-absorption measurements give a mean 
life of (1.15+0.14) x10™* sec. Debye temperatures 
of 400°K and 700°K were used for Li and LiF, 
respectively. The error given, about twice the 
statistical error of the data, includes estimates 
of the effect of other experimental uncertain- 
ties. 


EXPERIMENTS CONCERNING THE LOW-ENER- 
GY STATES OF THE O”* NUCLEUS. W. Zimmer- 
mann, Jr., Kellogg Radiation Laboratory, Cali- 
fornia Institute of Technology, Pasadena, Cali- 
fornia (Received December 1, 1958). 


Angular distributions have been measured for 
three groups of protons from the O"*(d,p)O”* re- 
action, those leaving O”* in its states at 0, 0.096, 
and 1.47 Mev. Deuteron energies of 1.74 and 2.50 
Mev in the laboratory system were used. The 
distributions of protons leaving O'* in its ground 
state and in its 1.47-Mev state are characteristic 
of stripping and indicate the formation of the 
ground state by an / =2 neutron and of the 1.47- 
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Mev state by an /=0 neutron. However, the dis- 
tribution of protons leaving O”* in its 0.096-Mev 
state does not lend itself to a stripping interpre- 
tation. 

It has been found that the )-decay of the 1.47- 
Mev state of O”*, following the formation of this 
state in the O'%(d, p)O"* reaction, proceeds 
mostly to the 0.096-Mev state. The mean life of 
the 0.096-Mev state has been measured by ob- 
serving the decay in flight of recoiling excited 
O”* nuclei and is found to be 1.75(140.16)x10~° 
second. These observations restrict the likely 
assignments of spin and parity for the 0.096-Mev 
state to 3/2* or 5/2°*. 


LEVELS IN Zr®: EXPERIMENTAL. S. Bjgrn- 
holm, O. B. Nielsen, and R. K. Sheline,* Institute 


for Theoretical Physics, University of Copenhagen, 


Copenhagen, Denmark (Received November 4, 
1957). 


The levels in Zr™ have been studied by analyzing 


the radiations of Nb® in magnetic and scintilla- 


tion spectrometers employing various coincidence 


techniques. Multipolarities of most of the transi- 
tions have been determined from internal con- 
version coefficients and K-L ratios. A decay 
scheme (I) for Nb® has been proposed which as- 
signs the following excited states in Zr™: 1752 
kev (0+), 2182 kev (2+), 2315 kev (5-), 3081 kev 
(4+), 3453 kev (6+), and 3595 kev (8+). Evidence 
has been discussed for a few weak additional 
transitions potentially involving three additional 
levels (decay scheme II). Following the sugges- 
tion of Ford, the levels in decay scheme I have 
all been interpreted as arising from the proton 
configurations (p,,)", (gy2)*, and (gy2 Py2). The 
half-life of the 3595-kev 8+ state has been ex- 
perimentally determined as 3x10~" sec, in good 
agreement with the half-life expected for a 141.5- 
kev E2 transition between 8+ and 6+ states, each 
involving a (g,.)* configuration. 

The relative population of the two 0+ states of 
Zr®™, both by de-excitation of the 2+ state of that 
nucleus and by the beta decay of Y®, indicates 
that these states result from highly mixed (p,,.)? 
and (g,.)* configurations. 

Hindrance factors for several transitions indi- 
cate that the other positive parity states are 
largely generated from the (g,,.)* configuration. 


"Now at Florida State University, Tallahassee, 
Florida. 
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LEVELS IN Zr™: THEORETICAL. B. F. Bayman, 
A. S. Reiner, and R. K. Sheline,* Institute for 
Theoretical Physics, University of Copenhagen, 
Copenhagen, Denmark (Received December 4, 
1958). 


An attempt is made to describe the seven levels 
of Zr®™ below 3.6 Mev in terms of the proton con- 
figurations (2py2)*, (2py2 1gy2), and (1gy2)*. The 
level positions and the compositions of the two 
0+ states are determined for Gaussian and Yukawa 
forces of various ranges and exchange characters, 
The experimental data are well reproduced for a 
reasonable choice of the force parameters, the 
best fit being obtained with a Serber exchange 
mixture and a range of about 1.5 fermis. The ex- 
perimental values of the lifetimes of the excited 
states can also be reconciled with these simple 
configurational assignments. The most serious 
discrepancy is in the half-life of the first excited 
(0+) state, which we calculate to be 1.35 x10" sec, 
as compared to the observed value of (6.0 + 1.5) 
x10-* sec. The remaining discrepancies in the 
energies and lifetimes are in the direction of the 
effects produced by a slight deformation of the 
Sr® core. 


— 
Present address: Chemistry Department, Florida 
State University, Tallahassee, Florida. 


TRITIUM PRODUCTION BY 6-Bev PROTONS. 
Lloyd A. Currie, Department of Chemistry, Penn- 
sylvania State University, University Park, Penn- 
sylvania (Received December 8, 1959). 


Tritium formation cross sections for 6-Bev 
protons on various targets between C and Pb are 
reported. No significant variations in tritium 
production cross sections are observed between 
2 Bev and 6 Bev. The problem of tritium pro- 
duction in nature is discussed briefly. 


PHOTOPRODUCTION OF K MESONS AND K- 
HYPERON PARITIES. Richard H. Capps, La- 
boratory of Nuclear Studies, Cornell University, 
Ithaca, New York (Received December 8, 1958). 


An investigation is made of the possibility of 
determining the intrinsic parities, relative to 
the nucleon, of K-A and K-= pairs by measuring 
the cross sections for photoproduction of these 
pairs from nucleons near threshold. The low- 
energy S- and P-wave contributions are calcu- 
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lated in a perturbation theory similar to that of 
Kawaguchi and Moravcsik. The question con- 
cerning the combination of nucleon and hyperon 
magnetic moments that is effective in the pro- 
cess is clarified. The calculated cross sections 
are quite different from those given in the long- 
photon-wavelength approximation. by Feld and 
Costa. For most reasonable choices of the 
hyperon anomalous magnetic moments, certain 
calculated ratios of the close-to-threshold cross 
sections corresponding to different K-A and K-= 
charge states are quite sensitive to the K-hy- 
peron parity. A model, in which the anomalous 
moments are the result of globally symmetric 
pion interactions, is used to illustrate the con- 
clusions. 


INNER BREMSSTRAHLUNG IN HYPERON DE- 
CAYS. S. Barshay* and R. E. Behrends, Brook- 
haven National Laboratory, Upton, New York 
(Received December 10, 1958). 


The angular distribution, energy spectrum, 
and total rate for the pion modes of hyperon de- 
cay with inner bremsstrahlung are calculated. 
The contributions from static anomalous mag- 
netic moments are also included. In principle, 
it is then possible to determine experimentally 
the one unknown parameter of the problem—the 
static anomalous magnetic moment of the hy- 
peron. 


*Present address: Institute for Theoretical Physics, 
Copenhagen, Denmark. 
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